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Abstract. We prove the immediate appearance of an explicit lower bound for 
continuous mild solutions to the full Boltzmann equation in the torus or a C 2 
convex domain with specular boundary conditions, under the sole assumption 
C*~) of regularity of the solution. We investigate a wide range of collision kernels, 

some satisfying Grad's cutoff assumption and others not. We show that this 
lower bound is exponential, independent of time and space with explicit constants 
depending only on the a priori bounds on the solution. In particular, this lower 
bound is Maxwellian in the case of cutoff collision kernels. A thorough study of 
characteristic trajectories, as well as a geometric approach of grazing collisions 
against the boundary arc derived. 
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1. Introduction 

This paper deals with the Boltzmann equation, which rules the behaviour of rar- 
efied gas particles moving in a domain U of lR d with velocities in M, d (d ^ 2) when 
the only interactions taken into account are binary collisions. More precisely, the 
Boltzmann equation describes the time evolution of f(t,x,v), the distribution of 
particles in position and velocity, starting from an initial distribution f Q (x,v) . 

We investigate the case where U is either the torus or a C 2 convex bounded domain 
and we assume the collisions against its boundary to be specular reflections. The 
Boltzmann equation reads 
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(1.1) Vt^O , V(x,v)eUxR d , d t f + vVJ = QU\f), 

\/(x,v)eUxR d , f(0,x,v) = f Q (x,v), 

with / being periodic in the case of hi = T d , the torus, or with / satisfying the 
specular reflections boundary condition if U is a C 2 convex bounded domain: 

(1.2) V(x, v) e dlA x R d , f(t, x, v) = f(t, x, TL x {v)). 

1Z X , for x on the boundary of U, stands for the specular reflection at that point of 
the boundary. One can compute, denoting by n(x) the outward normal at a point 
x on 8U, 

VveR N , K x (v) = v - 2(v ■ n(x))n(x). 

The quadratic operator Q(f, f) is local in time and space and is given by 

Q(/,/)= / B(\v-v*\,cos9)[f'fi-ff*]dv*da, 

where /', /*, fl and / are the values taken by / at v', v*, v'^ and v respectively. 
Define: 

/ = v+v* , \v-v»\ 

u 2 2 V V 



The collision kernel B ^ contains all the information about the interaction be- 
tween two particles and is determined by physics (see [2] or [3] for a formal derivation 
for the hard sphere model of particles). In this paper we shall only be interested 
with the case of B satisfying the following product form 

(1.3) B(|v- v,|,cos0) =$(|u-t;,|)&(cos0), 

which is a common assumption as it is more convenient and also covers a wide range 
of physical applications. Moreover, we shall assume that $ satisfies either 

(1.4) Vz e R, c$ \z\ y $(z) < C* \zf 

or a mollified assumption 

V|*|>leR, c$ N 7 $(z) < C$ |^| 7 
V|zj^leR, c$ ^ $(2) ^ C$ 

c$ and C$ being strictly positive constants and 7 in (-N, 1]. The collision kernel 
is said to be "'hard potential"' in the case of 7 > 0, "'soft potential"' if 7 < and 
"'Maxwellian'" if 7 = 0. 

Finally, we shall consider b to be a continuous function on 9 in (0,7r], strictly 
positive near 6 ~ ir/2, which satisfies 

(1.6) b (cos 6) sin d - 2 6 ~ 6 r (1+1/) 

0->O+ 

for 60 > and 1/ in (—00, 2). The case when b is locally integrable, v < 0, is referred 
to by the Grad's cutoff assumption (first introduce in [?]) and therefore B will be 
said to be a cutoff collision kernel. The case v ^ will be designated by non-cutoff 
collision kernel. 
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1.1. Motivations and comparison with previous results. The aim of this arti- 
cle is to show and to quantify the strict positivity of the solutions to the Boltzmann 
equation when the gas particles moves in a bounded domain. This issue has been 
tackled for a long time since it not only presents a great physical interest but also 
appears to be of significant importance for the mathematical study of the Boltzmann 
equation. 

Besides, our results only require some regularity on the solution and no further 
assumption on its local density, which was assumed to be uniformly bounded from 
below in previous studies (which is equivalent of assuming a priori either that there 
is no vacuum or that the solution is strictly positive). 

More precisely, we shall prove that continuous solutions to the Boltzmann equa- 
tion with angular cutoff in a C 2 convex bounded domain or the torus which have 
uniformly bounded energy satifies an immediate Maxwellian lower bound: 

Vt >0, 3p,9>0, Vt^t , W(x,v) eU xR d , f(t,x,v) > L/ 2 ^- 



The strict positivity of the solutions to the Boltzmann equation standing in the 
form of an exponential lower bound was already noticed by Carleman in [1] for the 
spatially homogeneous equation. In his article he proved that such a lower bound 
is created immediately in time in the case of hard potential kernels with cutoff in 
dimension 3. More precisely, the radially symmetric solutions he constructed in pQ 
satisfies an almost Maxwellian lower bound, 

Vt > t ,Vv e R 3 , f(t,v) > C ie - C2H2+E , 

GxiCi > for all to > and e > 0. His strategy was to start from an initial data 
uniformly bounded from below by a ball centred at the origin and then to use a 
spreading property of the collision operator to grow this ball bigger and bigger. 

The spreading strategy was followed by Pulvirenti and Wennberg in [12] to extend 
the latter inequality to solutions to the spatially homogeneous Boltzmann equation 
with hard potential and cutoff in dimension 3 with much more general initial data. 
The only assumption was a physically expected property of the solutions, that is the 
uniform boundedness of its density, its energy and its entropy. Their contribution 
was to get rid of the inital boundedness suggested in [1] by Carleman thanks to the 
use of a regularity property of the collision operator. This property allowed them to 
create immediately an "upheaval point" that they then spread with the method of 
Carleman. Besides, they obtain an exact Maxwellian lower bound of the form 

\/t^t ,Vv GM 3 , f(t,v) ^ C ie - C2lvl \ 

for all t > 0. 

Finally, the same approach than the one in [12] together with the use of the de- 
scription of solutions along characteristic trajectories has been used by Mouhot in 
[TP] to deal with the full Boltzmann equation with cutoff on the torus in d dimen- 
sion. He proved the immediate appearance of a Maxwellian lower bound, uniform 
in time and space, as long as the solutions has uniformly bounded density, energy 
and entropy (for the hard potential case) together with uniform bounds on higher 
moments (for the soft and Maxwellian potentials case). However, he also implicetly 
assumed that the initial data had to be bounded from below uniformly in space. 
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Mouhot also derived, in [10J, the same kind of results in the non-cutoff case in the 
torus. Under the same assumptions plus uniform boundedness assumption on higher 
Sobolev norms, he managed to prove the immediate appearance of an exponential 
lower bound of the form 

Vt ^ t ,V{x,v) ET d x R d , f(t,v) ^ C 1 (e)e- C2{£)lvlK+ \ 

for all to > 0, all e > and K = Kiv) with K(0) = 2 (thus recovering the cutoff 
case in the limit). In this case again it was implicitely assumed that the initial data 
is bounded from below uniformly in space. 

Our results extend those in [TU] by getting rid of any assumption on the local 
density and above all by dealing also with C 2 bounded convex domain. Furthermore, 
it does not assume any uniform boundedness on the initial data but requires the 
continuity of the solution to the Boltzmann equation. 

The quantification of the strict positivity, and above all the appearance of an ex- 
ponential lower bound has been seen to be of great mathematical interest thanks to 
the development of the entropy-entropy production method. This method (see [T5] . 
Chapter 3, and [IS]) provides a useful way of investigating the long-time behaviour 
of solutions to kinetic equations. Indeed, it has been successfully used to prove 
convergence to the equilibrium in non-pertubative cases for the Fokker-Planck equa- 
tion, [5j, and the full Boltzmann equation in the torus or in C 1 bounded connected 
domains with specular reflections, [6j. This entropy-entropy production method re- 
quires (see Theorem 2 in [6]) uniform boundedness on moments and Sobolev norms 
for the solutions to the Boltzmann equation but also an a priori exponential lower 
bound of the form 

f{t,x,v)^C ie - c ^\\ 

with q ^ 2. 

Therefore, the present paper allows us to prove that the latter a priori assumption 
is in fact satisfied for a lot of different cases (see [IP] . Section 5 for an overview). We 
also emphasize here that the assumption of continuity of the solution we have made 
does not shrink the range of applications since a lot more regularity is asked for the 
entropy-entropy production method so far. Moreover, our method, unlike the ones 
developed in [10] and [12], does not require uniform bounds on the local entropy or 
the local density of solutions, which neither is a requirement for the entropy-entropy 
production method (see [BJ, Theorem 2). 

To conclude we notice that our investigations require a deep and detailed under- 
standing of the geometry and properties of characteristic trajectories for the free 
transport equation. In particular, a geometric approach of grazing collisions against 
the boundary is derived. The existing strategies as well as our improvements are 
discussed in the next section. 

1.2. Our strategy. Our strategy to tackle this issue will follow the method intro- 
duced by Carleman, [JJ, together with the idea of Mouhot, [TU], to use Duhamel 
formula along characteristics. Roughly speaking we shall built characteristics in the 
torus or a C 2 bounded convex domain with specular reflections on the boundary and 
use the Duhamel formula to derive a lower bound for the solution to the Boltzmann 
equation. This lower bound should be first an "upheaval point" (as in [12] and [10] ) 
that we spread and grow bigger thanks to the spreading property of the collisional 
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operator (see Lemma 4.2). Finally, once the lower bound can be compared to an 
exponential one we reach the expected result. 

However, the existence of rebounds against the boundary leads to difficulties. We 
describe them below and point out how we shall overcome them. 

The first step of our study is to be able to construct the characteristic trajectories 
associated to the Boltzmann equation with specular reflections in a C 2 bounded 
convex domain. These trajectories are merely those of the free transport and so can 
be seen as the movement of a billiard ball inside the boundary of our domain. 

Such a free transport in a convex media has been studied in pE] (see also [13], [H] 
or [11] for geometrical properties) and has been used in kinetic theory by Guo, [8], 
or Hwang, [9J, for instance. Yet, the common feature in [3], [8] and [9] is that their 
assumptions on the boundary always lead to clear rebounds of the characteristic 
trajectories. That is to say, the absoption phenomenon of [1] or the electromagnetic 
field in |8j and [9] prevent the characteristics to roll on the boundary which is 
one of the possible behaviour we have to take into account in our general settings. 
As quickly mentionned in the introduction of [H], the behaviour at some specific 
boundary points is mathematically quite unexpected, even if that is of no physical 
relevance. We thus classify all the possible outcomes of a rebound against the 
boundary and study them carefully to analytically build the characteristics for the 
free transport equation in our domain U. 

Finally, we need to control the number of rebounds that can happen in a finite 
time. In [T3], Tabachnikov focuses on the footprints on the boundary of the tra- 
jectories of billiard balls and shows that the initial conditions leading to infinitely 
many rebounds on the boudary is a set of measure 0. We extend this to the whole 



trajectory (see Section 3.1, Proposition 3.4), not only its footprints on the boundary, 



allowing us to consider only finitely many rebounds in finite time and to have an 
analytic formula for the characteristics which we shall use throughout the article. 

Notice that all this study of the free transport equation will be done in the case 
of a merely C l bounded domain. 



Once the characteristics are understood, we follow the same idea as in [10J and 
focus on mild solutions to the Boltzmann equation along those trajectories (see 



Definition 2.2 and Definition 2.5). The Duhamel form of these solutions allows us 



to use the spreading property of the collision operator (see Lemma 4.2) which shall 
lead us to a global lower bound if we create an inital localised lower bound (this is 
the strategy developed in [12J and [TO]). 

Creating an "upheaval point" was achieved, in [12] and pi)], by using an iterated 
Duhamel formula and a regularity property of the collision operator relying of a 
uniform lower bound on the local density of the function. But the use of this 
property requires a uniform control along the characteristics on the density, the 
energy and the entropy of the solutions to the Boltzmann equation which is fine for 
the homogeneous case but made Mouhot consider initial datum bounded from below 
uniformly in space. Our way of dealing with the appearance of the "upheaval point" 
is rather different but include more general initial datum. We make the assumption 
of continuity of solutions to the Boltzmann equation and by compactness arguments 
we can construct a partition of our phase space where initial localised lower bounds 
exist, i.e., localised "upheaval points". 
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The case on the torus studied by Mouhot tells us that an exponential lower bound 
should arise immediately and therefore we expect the same to happen as long as 
the characteristic trajectory is a straight line. Unfortunately, the possibility for 
a trajectory to remain a line depends on the distance from the boundary of the 
starting point, which can be as short as one wants. This thought is the basis of our 
means for spreading the initial lower bound. We divided our trajectories into two 
categories, the ones which always stay close to the boundary (grazing collisions) and 
the others. For the latter we can spread our lower bound uniformly as noticed in 
[10J and for the case of grazing collisions we will thoroughly investigate the geometry 
of such trajectories, show that their velocity does not evolve a lot along time and 
mix it with the spreading property of the collision operator. Notice here that the 
convexity of U is needed for the study of grazing trajectories. 

The last behaviour to notice is the fact that specular reflections completely change 
velocities but preserve their norm. Therefore, the existence of rebounds against the 
boundary prevents us from using the spreading property of the collision operator in 
an inductive way as easily as in previous articles. However, spreading a lower bound 
only depending on the norm of the velocity can be done as specular reflections are 
norm-preserving. To conclude, our strategy is to spread the lower bound created 
at the "upheaval points" independently for grazing and non-grazing collisions up to 
the point when the lower bound we obtain depends only on the norm of the velocity. 
Mainly, our lower bounds will be balls in velocity that can be not centred at the 
origin and we shall grow them up to balls containing the origin. Then, the induction 
made in pU] on those centred balls will grow them up to a uniform exponential lower 
bound. 

Collision kernels satisfying a cutoff property as well as collision kernels with a 
non-cutoff property will be treated following the strategy described above. The only 
difference is the decomposition of the Boltzmann bilinear operator Q we consider in 
each case. In the case of a non-cutoff collision kernel, we shall divide it into a cutoff 
collision kernel and a remainder. The cutoff part will already be dealt with and 
a careful control of the L°°-norm of the remainder will give us the expected lower 
bound, smaller than a Maxwellian lower bound. 

1.3. Organisation of the paper. Section [2] is dedicated to the statement and the 
description of the main results proved in this article. It contains four different parts 



Section 2.1 defines all the notations which will be used throughout the article. 



As mentionned above, we shall investigate in details the characteristics and the free 



transport equation in a C 1 bounded domain. Section 22 mathematically formulate 
the intuitive ideas of trajectories. 



The last subsections, 2.3 and 2.4, are dedicated to a mathematical formulation 



of the results related to the lower bound in, respectively, the cutoff case and the 



non-cutoff case, described in Section 1.2 It also defines the concept of mild solutions 



to the Boltzmann equation in each case. 



As said when we described our strategy (Section 1.2), we are going to study the 
free transport equation and the different important properties of the characteris- 
tics. Section [3] formulates these issues, investigates all the different behaviours of 



rebounds against the boundary (Section |3T), builds the characteristics and derives 



their properties (Section [3^2) and solves the free transport equation (Section 3.3). 
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Section [4] focuses on the proof of Theorem 2^3 about the Maxwellian lower bound 
in the cutoff case. It is divided into three subsections. 

Following our strategy, Section |4.1| creates the localised "upheaval points" and 



Section 4.2 proves the immediate appearance of a lower bound only depending on 
the norm of the velocity ( Proposition 2.4). We focus on the case of trajectories that 
are not entirely grazing ones (Section 



and finally we proved Proposition 2.4 



2.3) 



4.2.1), the case of grazing ones (Section 4.2.2) 



(Section 4.2.3) 



Section 4.3 derives the immediate Maxwellian lower bound (proof of Theorem 



Finally, we deal with non-cutoff collision kernels in Section [5] where we prove the 



immediate appearance of an exponential lower bound (Theorem 2.6). The proof 
follows exactly the same steps than the one in the case of cutoff collision kernel and 



so is divided into Section 5.1, where we construct a lower bound only depending 



on the norm of the velocity, and Section 5^2, where we derive the exponential lower 
bound. 

2. Main results 
We begin with the notations we shall use all along the article 



2.1. Notations. First of all, we denote (■) = yl + |-| 2 and y + = max{0,y}, the 
positive part of y. 

This study will hold in specific functional spaces regarding the v variable that we 
describe here and use throughout the sequel. Most of them are based on natural 
Lebesgue spaces L p v = LP {R N ) with a weight: 

• for p G [1, oo] and q G R, LP qv is the Lebesgue space with the following norm 

u/ik. = \\f(v) q \y v , 

• for p G [1, oo] and k G N we use the Sobolev spaces W^ ,p by the norm 

1 Vp 



k:2 



with the usual convention = W, 



In the sequel of this study, we are going to need bounds on some physical observ- 



ables of solution to the Boltzmann equation (1.1). 



We consider here a function f(t, x, v) ^ defined on [0, T) xU x R. N and we recall 
the definitions of its local hydrodynamical quantities. 



its local energy 

e f (t,x) 
its local weighted energy 

e' f (t,x) 
where 7 = (2 + 7) + , 



\v | 2 f(t, x, v)dv, 



\v | 7 fit, x, v)dv, 
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its local LP norm (p G [1, +00)) 

l p f (t,x) = \\f(t,x,.)\\ L „ 

its local W 2,oa norm 

Wfit.x) = \\f(t,x,-)\\ w > 



2,00 . 



Our results depends on uniform bounds on those quantities and therefore, to 
shorten calculations we will use the following 

Ef — sup ef(t,x) , E'f — sup e'f(t,x), 
(t,x)e[o,T)xU ' (t,x)e[o,T)xU 

L P f = sup E(t,x) , Wf — sup Wf(t,x). 

(t,x)<=[0,T)xU (t,x)e[0,T)xU 



In our theorems we are giving a priori lower bound results for solutions to (1.1) 
satisfying some properties about their local hydrodynamical quantities. Those prop- 
erties will differ depending on which case of collision kernel we are considering. We 
will take them as assumptions in our proofs and they are the following. 

• In the case of hard or Maxwellian potentials with cutoff (7^0 and v < 0): 

(2.1) E f < +00. 

• In the case of a singularity of the kinetic collision kernel (7 G (— N, 0)) we 
shall make the additional assumption 

(2.2) L p f < +00, 

where p 7 > N/ (N + 7). 

• In the case of a singularity of the angular collision kernel [y G [0, 2)) we shall 
make the additional assumption 

(2.3) W f < +00, E'f < +00. 

As noticed in [TO], in some cases several assumptions might be redundant. 

Besides, in the case of the torus with periodic conditions or the case of bounded 



domain with specular boundary reflections , solutions to (1.1), also satisfy the fol- 
lowing conservation laws (see [2] , [3] or [15] for instance) for the total mass and the 
total energy: 



(2.4) 3M, E > 0, Vt G 



Su fn d /(*» x ' v "> dxdv = M ' 
lu JvL d M 2 x , v ) dxdv = E. 



2.2. Results about the free transport equation. Our investigations start with 
the study of the characteristics of the free transport equation. We only focus on the 
case where U is not the torus (the characteristics in the torus being merely straight 
lines) but we will use the same notations in both cases. This is achieve by the 
following theorem. 
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Theorem 2.1. LetU be an open, bounded and C 1 domain in R . 
Letu :U xR N — >R be C 1 in x and in Lz „. 

Then the free transport equation with specular reflections, where TZ X stands for the 
specular reflection at a point x on the boundary, which reads 

(2.5) Vt^O , V(x,v) G U x R N , d t u(t,x,v) + D x (v)(u)(t,x,v) = 0, 

(2.6) V(i,!))GWxl N , u(0,x,v) =u (x,v), 

(2.7) V(i,i))eA + , u(t, x, v) = u(t, x, 1Z x (v)), 

has a unique solution u : l + xWxf w — > R which is C 1 in time, admits a directional 
derivative in space in the direction of v and is in L\ v . 

Moreover, for all (t,x,v) in R + x U x R N , it exists Xfi n (t,x,v), Vfi n (t,x,v) and 
tf in (t,x,v) (see Definition 3.6) such that 

u(t,X,v) = U (Xfin ~ (t - t fin )v fin ,V fin ) . 



2.3. Maxwellian lower bound for cutoff collision kernels. The final theorem 
we prove in the case of cutoff collision kernel is the immediate appearance of a 
uniform Maxwellian lower bound. We use, in that case, the Grad's splitting for the 
bilinear operator Q such that the Boltzmann equation reads 

Q(g,h) = / $ (\v - v*\) b (cos#) [h'gl - kg*} dv*da 

= Q+{g,h)-Q-{g,h), 
where we used the following definitions 

Q + (g,h) = / $(\v-v*\)b(coa9)h'£dv*d(T, 

(2.8) Q~(g,h) = n b {$*g{v))h = L[g](v)h. 



In Section[3jwe prove, Proposition (3.8), that we are able to construct the charac- 
teristics (X t (x, v), Vt(x, v )), for all (t, x, v) in IR + xWx R d , of the transport equation. 
Thanks to this Proposition we can define a mild solution of the Boltzmann equation 
in the cutoff case. This weaker form of solutions is actually the key point for our 
result and also gives a more general statement. 



Definition 2.2. Let f Q be a measurable function, non-negative almost everywhere 
onU x R d . _ 

A measurable function f = f(t,x,v) on [0,T) x U x R d is a mild solution of the 
Boltzmann equation associated to the initial datum f (x,v) if 

(1) f is non-negative onU x R d , 

(2) for every (x,v) inU x R d : 

t .— ► L[f(t, X t {x, v), .)](V t (x, v)),t^ Q + [f(t, X t {x, v), ■), fit, X t {x, v), -)]{V t {x, v)) 

are in L} oc ([0,T)), 

(3) and for each t G [0, T), for all x and v G R d , 
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- / L[f(s,X s (x,v),-)](V s (x,v))ds 





(2.9) f(t,X t (x,v), V t (x,v)) = f (x,v)exp 

+ J exp (- L[f( y s',X s/ (x,v),-)](V s/ (x,v)) ds' 

Q + [f{s,X s (x,v),-),f( S ,X s {x,v),-)}(V s (x,v))ds. 



Now we state our result. 



Theorem 2.3. LetlA be T d or a C 2 open convex bounded domain in M. d and let Jq be 
a continuous function onlA x M. d . Let B = Qb be a collision kernel satisfying (1.3), 



with $ satisfying (1.4) or (1.5) and b satisfying (1.6) with u < 0. Let f(t,x,v) be a 
mild solution of the Boltzmann equation in U x W 1 on some time intervalle [0, T), 
T G (0, +oo], which satisfies 

• f is continuous on [0, T) x U x Mr and f(0,x,v) = fo(x,v); 



if $ satisfies (1.4) with 7^0 or if Q satisfies (1.5), £/ien / satisfies (2.1); 



z/$ satisfies (1.4) 7 < 0, i/ien / satisfies (2.1) and (2.2) 



Taen /or all r G (0,T) i/iere exists p > and 6* > 0, depending on r, Ef (and Vf 
if $ satisfies (1.4) u;zi/i 7 < Oj, stxc/i i/iat /or a// t G [r, T) i/ie solution f is bounded 



from below, almost everywhere, by a global Maxwellian distribution with density p 
and temperature 6, i.e. 



Vi G [r,T), V(x,v) eUx 



f(t,x,v) ^ 



P 



{2ti6) n / 2 



As stated in Subsection 1.2, the main result to reach Theorem 2.3 is the construc- 



tion of an immediate lower bound only depending on the norm of the velocity: 



Proposition 2.4. Let f be the mild solution of the Boltzmann equation described 
in Theorem \2.3\ 

For all < r < T there exists ry, oo(t) > such that 



Vi G [r/2,r], W(x,v) G W x 



f{t,X,v) > Oo(r)lB{0,r v )( v ) 



ry and ao(r) on/y depending on r, Ef (and LF^ 1 if $ satisfies (1.4) wraia 7 < 0j. 



2.4. Exponential lower bound for non-cutoff collision kernels. In the case 



of non-cutoff collision kernels (0 ^ ^ < 2 in (1.6)), Grad's splitting does not make 



sense anymore and so we have to find a new way to define mild solutions to the 



Boltzmann equation (1.1). The splitting we are going to use is a standard one and 
it reads 



Q(g,h) 



$ (\v - f*|) b 

K iV xS JV-l 

Ql(g,h)-Qt(g,h), 



[cos9) [h'gl — hg*} dv*da 
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where we used the following definitions 

Ql(g,h) = [ <S>(\v-v*\)b{cos6)g',{ti -h) dv*da, 

(2.10) Q 2 b (g,h) = -(( <S>(\v-v*\)b{cos6)[gi-g*} dv*do)h 

\Jrn x§ n-i J 

= S[g](v)h. 

We would like to use the properties we derived from the study of collision kernels 
with cutoff. Therefore we are going to go further into the splitting we will use for 
the bilinear operator Q. 

For e in (0, 7r/4) we define a cutoff angular collision kernel 



6^(cos0) = b{cos9)l m>£ 

and a non-cutoff one 



6f co (cos#) = b (cos#) 1 



Considering the two collision kernels B^° = $6^° and B^ co = §b^ co , we can 



combine Grad's splitting (2.8) applied to with the non-cutoff splitting (2.10) 
applied to B^ co . This yields the splitting we shall use to deal with non-cutoff 
collision kernels, 



(2.11) Q = Q+_Q- + Ql-Q 



where we use the shortened notations Qf = Qf C o an d Ql^ 2 = Q 1 Jnco- 



Thanks to the splitting (2.11) and the study of characteristics announced in Sec- 
tion 



2.2 



we are able to define mild solutions to the Boltzmann equation with non- 
cutoff collision kernels. This is obtained by considering the Duhamel formula asso- 
ciated to the splitting (2.11) along the characteristics (as in the cutoff case). 



Definition 2.5. Let /o be a measurable function, non-negative almost everywhere 
onU x R d . _ 

A measurable function f = f(t,x,v) on [0,T) x U x W 1 is a mild solution of the 
Boltzmann equation with non- cutoff angular collision kernel associated to the initial 
datum fo(x, v) if there exists < e < ir/4 such that for all < e < e : 

(1) / is non-negative onU x M. d , 

(2) for every (x,v) inU x M. d : 

t ^ L e [f(t, X t (x, v), •)] (V t (x, v)),t^ Q+[f(t, X t (x, v),-), f(t, X t (x, v), -)}(V t (x, v)) 

t — ► S £ [f{t, X t (x, v), ■)] (V t (x, v)),t^ Ql[f(t, X t (x, v), ■), f(t, X t (x, v), -)}{V t (x, v)) 

are in L} oc ([0,T)), 

(3) and for each t G [0, T), for all x and v G M. d , 



12 



M. BRIANT 



(2.12) 

f(t,X t (x,v),V t (x,v)) = f (x,v)exp 



- / {L £ + S £ )[f(s,X s {x,v),-)](V s (x,v)) ds 



+ / exp 
'o 



(L £ + S e ) [f(s',X s ,(x,v),-)}(V s ,(x,v)) ds' 
(Qt + Ql) [f(s, X s (x, v),-), f{s, X s (x, v), -)}(V s (x, v)) ds. 



Now we state our result. 



Theorem 2.6. Let U be T d or a C 2 open convex bounded domain in M, d and fo be 



a continuous function onlA x M. . Let B = $6 be a collision kernel satisfying (1.3), 



with $ satisfying (1.4) or (1.5) andb satisfying (1.6) with v in [0,2). Let f(t,x,v) be 
a mild solution of the Boltzmann equation inU xtf on some time intervalle [0,T), 
T G (0, +oo], which satisfies 

• f is continuous on [0, T) xWxtf and f(0,x,v) = fo(x,v); 



if $ satisfies (1.4) with 7 ^ or if $ satisfies (1.5), then f satisfies (2.1) 



and (2.3); 

• if& satisfies ( Jl~4| ) urat/i 7 < 0, then f satisfies ( |2~Tj ), ( |2~2"| and ( |2~3| ) 

Taen for all r G (0, T) and /or any exponent K such that 



K > 2 



(2 + a?) 

log2 



i/iere exists C\,C2 > 0, depending on r, 



(1.4) u>zta 7 < Oj, sitc/i taat 

Vt G [r, T), V(z, v)eUx R d , f(t, x, v) ^ de^M* . 
Moreover, in the case v = 0, one can take K = 2 (Maxwellian lower bound). 



We emphasize here that, in the same spirit than in the cutoff case, the main part 



of the proof will rely on the establishment of an equivalent to Proposition 2.4 for 
non-cutoff collision kernels. 



Remark 2.7. Throughout the paper, we are going to deal with the case where U is a 
C 2 convex bounded domain since it is the case where the most important difficulties 
arise. However, ifU = T d , we can follow the same proofs by letting the first time of 
collision with the boundary to be +00 (see Section^ and by making the definition 
that the distance to the boundary (which does not exist) is +00 (which rules out the 
case of grazing trajectories). 
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3. The free transport equation: proof of Theorem |2~T 



In this section, we study the transport equation with a given initial data and 
boundary condition in a bounded domain U. We will only consider the case of 
purely specular reflections on the boundary dtl. Those kind of interaction cannot 
occur for all velocities at the boundary. Indeed, for a particle to bounce back at the 
boundary, we need its velocity to come from inside the domain U. To express this 
fact mathematically, we define 

A+ = {( x ,v) G dU x R N : v ■ n(x) ^ 0} , 

where we denote by n(x) the exterior normal to dU at x. 

Consider u : U x R N — > R which is C 1 in x G U and L 2 (U x R N ) = L 2 XV . We 



are interested in the problem stated in Theorem 2.1, (2.5) — (2.7). 

If D x (v)(u) denotes the directional derivates of u in x in the direction of v we 
have, in the case of functions that are C l in x, 

D x (v){u) = v ■ V x u. 

Therefore, instead of imposing that the solution to the transport equation should be 
C 1 in x, we reformulate the problem with directional derivatives, and that allowed 
us to prove the following theorem. 

Physically, the free transport equation means that a particle evolves freely in U at 
a velocity v until it reaches the boundary. Then it bounces back and moves straight 
until it reaches the boundary for the second time and so on so forth up to time 
t. The method of characteristics is therefore the best way to link u(t,x,v) to u 
by just following the path used by the particle, backwards fromt to (see Figjl]). 
This method has been used in jH] on the half-line and in [1], [0], for instance, in 




t t (t v ») 



FIGURE 1 . Backward trajectory with standard rebounds 



the case of convex media. However, in both articles they only deal with finite, or 
countably many, numbers of rebounds in finite time. Indeed, the electrical field 
in [8] and [9] makes the particles always reach the boundary with v ■ n(x) > 
and [1] has a specular boundary problem with an absorption coefficient a G [0, 1): 
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u(t, x,v) = au(t,x,lZ x (v)). Therefore, in the case the particle arrives tangentially 
to the boundary, i.e. v ■ n(x) = 0, we have lZ x (v) = v and so u(t,x,v) = 0. This 
vanishing property allowed the authors not to care about special cases where the 
particle starts to roll on the boundary since the value of the function is given by the 
boundary conditions. 

Another way of looking at the characteristics method is to study the footprints 
of the trajectories on the boundary. This problem, as well as the possibility of 
having infinitely many rebounds in a finite time, has been tackled by Tabachnikov 
in [13] . Tabachnikov only focused on boundary points since the description of the 
trajectories by only considering their collisions with the boundary raises a symplectic 
property and a volume-preserving transformation. Such properties allowed him to 
show that the set of points on the boundary that leads to infinitely many rebounds 
in finite time is of measure ([13], Lemma 1.7,1). Unfortunately, in our case we 
would like to follow the characteristics and the study of trajectories only via their 
footprints on the boundary is no longer a volume-preserving transformation. 

In our case we do need to follow the path of a particle along the characteristics of 
the equation to know the value of our function at each step. If the particle starts to 
roll on the boundary (see Fig|2| we require to know for how long it will do so. The 
major issue is the fact that v ■ n(x) = do not tell us much about the geometry of 
dU at x and the possibility, or not, for the particle to keep moving tangentially to 
the boundary. Moreover, some cases lead to non physical behaviour since the sole 
specular collision condition implies that some pairs (x,v) G dU x M, N can only be 
starting points, they cannot be generated by any trajectories (see Figj3]). This case 
is mentioned quickly in the first chapter of [H] but not dealt with. 




Therefore, in order to prove the well-posedness of the transport equation (2.5) — 
(2.7), we follow the ideas developed in [8] and [9], which consist in studying the 
backward trajectories that can lead to a point (t,x,v), combined with the idea of 
countably many collisions in finite time used in [1]. However, we have to deal with 
the issues described above and to do so we introduce a new classification of possible 
interactions with the boundary (see Definition 3.1). We also manage to extend the 
result of [13], in terms of pair (x,v) leading to infinitely many rebounds in finite 
time, to the whole domain U (Proposition 3.4). To do so we link up the study on 
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FIGURE 3. Backward trajectory that reaches an end 

the boundary made in [13] with the Lebesgue measure on U by artificially creating 
volume on dU thanks to time and a foliation of the domain by parallel trajectories. 

The section is divided as followed. First of all we shall describe and classify the 
collisions with the boundary in order to describe very accurately the backward trajec- 
tories of a point (x, v) in dli x R . We will name trajectory or characteristic any so- 
lution (X(t, x, v), V(t, x, v)) satisfying the initial condition (X(0, x, v), V(0, x, v)) = 




(x,v), the boundary condition ( |2.7[ ) and satisfying, in U, 

= V 
= 0. 

This will give us an explicit form for the characteristics and so allow us to link 
u(t, x, v) with uo(x*, v*), for some x* and v*. Finally, we shall show that the function 
we constructed is, indeed, a solution to the transport equation with initial data uq 
and specular boundary condition and that such a solution is unique. 

3.1. Study of rebounds on the boundary. As mentionned in the introduction 
of this section, when a particle reaches a point at a boundary with a velocity v it can 
bounce back, FigjTJ keep moving straight, FigfJJor, and this is physically unexpected, 
stop moving because the specular reflection does not allow it to do anything else, 
Fig|3j The next definition gives a partition of the points at the boundary which 
takes into account those properties. 

Definition 3.1. We define here a partition of dU x M. N that focuses on the outcome 
of a collision in each of the sets. 

• The set coming from a rebound without rolling 

^rebounds = {(,X,v) G OU X M. N : V • n{x) < 0} . 

• The set coming from rolling on the boundary 

^roiling = {(x, v) e dU x R N : v u(x) = and 35 > 0, Vt G [0, 5), x-vteU}. 

• The set of only starting points 

tt stop = {(x,v) e dU x R N : vn(x) = and V5 > 0, 3t e [0,5], x-vt^U}. 
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• The set coming from straight line 

ttu n e = {(x,v) G dU x R N : v n(x) > 0} . 

One has to notice that any point of flu ne indeed comes from a straight line arriving 
at x with direction v since U is open and is C 1 (so there is no cusp). 
In order to understand the behaviour expected at VL stop we have the following propo- 
sition. The proof of it really shows what are specular boundary reflections. 



Proposition 3.2. If we have (x,v) in fl s top then there is no trajectory with specular 
boundary reflections that leads to (x, v) . 



Proof of Proposition 3J2 , Let us assume the contrary, that is to say (x, v) is in VL stop 
comes from a trajectory with specular boundary reflection. 

We have that (x, v) belongs to dtl x M. N and so if (x, v) comes from a straight 
line it can only be (by definition of trajectories) a line containing x with direction v 
which means that (x, v) comes from {(x — vt,v), t G [0, T]}, for some T > 0. But a 
trajectory is necessarily in U and this is in contradiction with the definition of ^stop- 
Therefore, (x, v) must come from a rebound after a straight line trajectory. But 
again we obtain a contradiction because the velocity before the rebounds is lZ x (v) = 
v and the backward trajectory is the one studied above. □ 



Now we have our partition of points on the boundary of U, we are able to generate 
the backward trajectory associated to a starting point (x, v) in tl x M. N . The first 
step towards its resolution is to find the first point of real collision (if it exists) that 
generates (x,v) (see Fig{TJ. The next proposition-definition proves mathematically 
what the figure shows. 



Proposition 3.3. Let U be an open, bounded and C 1 domain in K . Let (x,v) be 
inlA x M. N , then we can define 

t min (x, v) = max {t ^ : x — v s G U, V ^ s ^ i j . 

Moreover we have the following properties: 

(1) for allt in (0, t min (x, v )), ifx — vt hits dU then (x — vt,v) belongs to Q r oiiing- 

(2) t min (x,v) = if and only if (x,v) belongs to Q stop U ^rebounds- 

(3) x - vt min (x,v) belongs to VL stop U VL rebounds . 



Property 1. emphasises the fact that if, on the straight line between x and x — 
vt m i n (x,v), the particle hits the boundary it will not be reflected and so just rolls 
on. Then property 2. tells us than t min ( x ^ is always strictly positive except if (x, v) 
does not come from any trajectory of a particle or if it is the outcome of a rebound 
without rolling. Finally, property 3. finishes the study since at x — vt min (x,v) the 
particles either come from a reflection (case fl re bounds), and we can keep tracking 
backwards, or started its trajectory at x — vt min (x,v) (case fl s to P )- 
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Proop of Proposition |3.3[ First of all we have that IA is bounded and so there exists 
R such that U C B(0, R), the ball of radius R in R N . 
Then we notice that belongs to 

A(x, v) = {t ^ : x - vs E U, V < s < t) . 

Therefore A(x,t>) is not empty. Moreover, this set is bounded above by 2i?/ ||i> || 
since for all t in -A (a;, t>) 

R> \\x — vt\\ ^ t \\v\\ — \\x\\ . 

Therefore we can talk about the supremum t m i n (x,v) of A(x,v) and there exists an 
increasing sequence (t n ) n( z N in A(x,v) that tends to t min (x,v). As U is closed we 
have that x — vt min (x,v) belongs to U. Then, if ^ s < t min (x,v) there exists n 
such that ^ s ^ t n and so, by the property of t n , x — v s is in U. This conclude 
the fact that t min (x, v) belongs to A(x, v) and so is a maximum. 

No we turn onto the proof of properties. 

Let (x, v) be in U and < t < t min (x, v) such that x — vt belongs to dlA. Then for 
all < t\ < t < t 2 < t min (x, v), x — vti and x — vt 2 are in U and so, by the definition 
of an exterior normal to a surface we have 

[(x — vt) — (x — vti)] ■ n(x — vt) ^ and [(x — vt) — (x — vt^)] ■ n(x — vt) ^ 0, 

which gives v ■ n(x — vt) = 0. 

Moreover, since t 2 belongs to A(x,v), for all s in [0,t 2 — t], (x — vt) — vs is in U, 
which means that (x — vt, v) belongs to ^rolling- 

Property 2. is direct since if t min (x, v) — then for all t > 0, there exists < s ^ t 
such that x — vs does not belong to U and then w • n(x) ^ 0. So (x,v) belongs to 
^rebounds, if v • n(x) > 0, or in fi stop . 

Finally, property 3. is straightforward since x — vt min (x, v) is in dlA (because U is 
open) and since for all ^ t ^ t min (x, v), x—vt is in U. thus [{x — vt min ) — (x — vt)}- 
n(x — vt min (x, v)) ^ 0, which raises v ■ n(x — vt min (x, v)) ^ O.Then, by the property 
of A(x,v) and the fact that t m i n (x, v) is its maximum, we have that either (x — 
vt min (x, v),v) belongs to VL rebounds or belongs to Q stop . □ 



Up to now we have been focusing only on the case of the first possible collision 
with the boundary. In order to conclude the study of rebounds for any given char- 
acteristics we have to, in some sense, count the number of rebounds without rolling 
that can happen in finite time. This is the purpose of the next proposition. 



Proposition 3.4. Let U be a C 1 open, bounded domain in M, N and let (x, v) be in 
U x R N . 

Then for all t ^ the trajectory finishing at (x,v) after a time t has at most a 
countable number of rebound without rolling. 

Moreover, this number is finite almost surely regarding the Lebesgue measure on 
UxR N 



Proof of Proposition 3.4, The fact that there is countably many rebounds without 
rolling comes directly from the fact that t m i n (x,v) > except if (x,v) is a starting 
point (and then did not move from to t) or if (x, v) is the outcome of a rebound (and 
so comes from (x,lZ x (v)) which belongs to Qu ne , implying that t m i n (x,7l x (v)) > 0). 
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Therefore the time between two distinct rebounds without rolling is always strictly 
positive which gives us that the number of rebounds without rolling is at most 
countable. 

Now we shall prove that the set of points in W x R N which lead to an infinite 
number of rebounds in a finite time is of measure 0. To do so, we first need some 
definitions. The measure \x in 14 x M. N is the one induced by the Lebesgue measure 
and we denote by A the measure on dlA x M. N (see section 1.7 of [T3]). 
We will also denote 

VI = {(x, v) G 14 x (R — {0}) coming an infinite number of rebounds} , 

— {{ x i v ) 914 x (R — {0}) coming from an infinite number of rebounds} . 

We know, [T3] Lemma 1.7.1, that X(flg) = and we are going to establish a link 
between the measure of fl and the one of fig. Those two sets do not live in the same 
topology and so we build a function that artificially recreates the same topology via 
time. 

Because 14 is bounded we can find a point X outside 14 and a time TV > 
such that for all x in U and v in WL N — {0}, (x — Tmv/ \\v\\) does not belong to U. 
Furthermore, in the same way as for t m i n (x, v), we can define, for (x, v) in 14 x R N , 

T(x, v) = min{t ^ : x + vt G dK\. 

Then we have all the tools to build up a link between Q and Qg. Let us define the 
following function which is clearly C 1 . 

F : [0, T M ] x R n x (R n - {0}) — ► R N x (R* - {0}) 

(t,x,v) .— > + 

There is a set where F has good properties: 

B = {{t,x,v) : x G v e (R N — {0}), t G [0,T(z,u))} . 

F is injective on the set B. 
Indeed, if (t,x,v) and (t*,x*,v*) are in B such that F(t,x,v) = F(t*,x*,v*) then 
v = v* and x + tv / ||u|| = x* + ||u||. 
Let assume that t* > t, therefore we have that 

x = x * + (t* _t) e aw 

1 1 17 1 1 

and thus t* — t ^ T(x*,v). However, t* ^ T(x*,v) so we reach a contradiction and 
t* ^ t. By symmetry we have t = t* and then x = x*. 

Finally we have that fl = F (B D ([0,7V] x ^a))- Indeed, if (t,x, u) belongs to 
Bfl([0, Tm] x fig) then F(t, x, u) = (x + tv / \\v\\ , u) and x + tf / ||f || is in U and its 
first rebound backward in time is (x, v) which lead to infinitely many rebounds in 
finite time. Therefore 

x + tj-j G 

The converse is direct, by considering the first collision with the boundary of the 
backward trajectory starting at (x, v) in fl. 

All those properties allow us to compute /i(fi) by a change of variable in B fl fig. 
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//(fi) = fi(F(Bn{[Q,T M ]xn d ))) 

^F(Bn([o,T M ]xn a )){x,v)dxdvdt 

UxR N 

|Jac(F _1 )| d\(x,v)dt 



* Bn([o,T M ]xn a ) 
< T M sup (jJac^ 1 )!) X(n d ) = 0. 

[0,T M ]xU 

□ 



3.2. Description of characteristics. In the previous section we derived all the 
relevant properties of when, where and how a trajectory can bounce against the 
boundary of hi. As said above, the characteristic starting from a point (t,x,v) in 
1R + x U x M. N is the backward trajectory satisfying specular boundary reflections 
that leads to (x, v) in time t. Basically, it consists in a straight line as long as it 
stays inside U or it rolls on the boundary. Then it reaches a boundary point where 
it does not move any more (fl s t op ) or bounces back (fl re bounds) ■ 

Thanks to Proposition 3.4 we can generate the countable (and almost surely finite) 
sequence of collisions with the boundary associated to the future point (x,v). We 
shall construct it by induction. We consider (x,v) in W x K . 

• Step 1: initialisation: we define 



x (x,v) 
v (x,v) 
t (x,v) 



x, 
0. 



Step 2: induction: if (xk(x,v),Vk(x,v)) € fi s t op then we define 



x k+ i{x,v) 
v k+1 (x,v) 
t k+1 (x,v) 



x h (x,v), 
v k (x,v), 
+oo, 



if (xk(x,v),Vk(x,v)) ^ fl s top then we define 

x k+ x(x,v) = x k {x,v) - v k (x,v)t min (x k (x,v),v k (x,v)), 
v k+1 (x,v) =K Xk+l{x , u) (v k (x,v)), 
t k+1 (x,v) = t k (x,v) +t min (x k (x,v),v k (x,v)). 

Remark 3.5. Let us make a few comments on the accuracy of the sequence we just 
built. 



(1) Looking at Proposition 3.3 , we know that at each step (apart from 0) we nec- 
essary have that (x k (x,v),v k (x,v)) belongs to either VL stop or ^rebounds an d 
so the characteristic stops for ever (case 1 in induction) or bounces without 
rolling and start another straight line (case 2). Thus the sequence of foot- 
prints defined above captures the trajectories as long as there are rebounds 
and then becomes constant once the trajectory reach a stopping point. 
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(2) If t m i n (xk(x,v),Vk(x,v)) = for some k > then, by properties 2. and 3. 



of Proposition 3.3, we must have (x k (x,v),v k (x,v)) G Vt stop (since v k (x,v) is 
the specular reflection atx k (x,v) of Vk-i(x,v) and (xk(x,v),Vk-i(x,v)) is in 
fl re bounds UQ stop ). Thus, (t k (x, v ))& e N strictly increasing as long as it does 
not reach the value +oo, and then it remains constant. 

Finally, it remains to connect the time to those quantities. In fact, the time will 
determine how many rebounds can lead to (x, v) in a time t. The reader must 
remember that the backward trajectory can lead to a point in fl s t op before time t. 

Because the characteristics method can help us knowing the value of the solution 
to the transport equation at a given point thanks to its value at any other point of 
the trajectory, the next definition links a triplet (t, x,v) to the first rebound of the 
trajectory that leads to (x, v) in a time t. 



Definition 3.6. LetlA be an open, bounded and C 1 domain in K . 
Let (t,x,v) be in IR + x U x M. N . Then we can define 

max{fc G N : tk(x,v) ^ t}, if it exists, 
+oo, if (tk(x,v))k is bounded by t. 

The last rebound is then define by 

• ifn(t,x,v) < +oo and tnfcx^+i = +oo, then 

X fin (t) X , V) X n (f t x t v) \X , Vj , 
Vf in (t,X,v) = V n (f iX)V )(x,v), 

t fin (t,x,v) = t, 

• ifn(t,x,v) < +oo and t n (f )XiV -)+i < +oo, then 

Xfin(t)X,V) X n (f t x t v) \X , Vj , 
v fin (t,x,v) = V n ( t>x>v )(x,v), 
tfin\f , iX,V) tn(t,x,v) \X, V) , 

• ifn(t,x,v) = +oo, then 

{Xf in (t,x,v) = lim Xk(x,v), 
k— >+oo 
Vf in (t,x,v) = lim v k {x,v), 
k— 5>+oo 
tf in (t,x,v) = lim t k (x,v). 
k— S>+oo 

Remark 3.7. Let us make a few comments on the definition above and the existence 
of limits. 

(1) After the last rebound, occuring at t n (t, x ,v)> the backward trajectory can only 
be a straight line during the time period t — t n (t, x> v) ( see That is 
why we defined tf in (t,x,v) = t n (t,x,v) if we reached a point on Vt r ebounds and 
tfi n (t,x,v) = t if the last rebounds reaches Q s top (the trajectory can only start 
from there). 

(2) In the last case of the definition, we remind the reader that (tk{X)V))keti ^ s 
strictly increasing and so converges if bounded byt. But then, t min (xk{x, v),Vk(x, 
tends to as k goes to +oo. Moreover, (xk(x,v)) k( z N is in the compact dU 
and so has a convergent subsequence. Then because (||i>fc(x, v ) |))fceN ^ s con- 
stant and x k (x, v) = x k -i(x, v) — t min (xk(x, v),V)~(x, v))vk(x, v), we have that 
(xk{x, v))fc S N converges and thus (v k (x,v))keN as well. 



n(t, x, v) 
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(3) The last case in Definition 3.6 almost surely never happens, has proved in 



Proposition 3.4 



To conclude this study of the characteristics we just have to make one more 
comment. We studied the characteristics backward in time because it is easier in 
order to construct a solution to the free transport equation. However, it is easy 
to prove (just require the inductive construction of Vk and Xk) that the forward 
trajectory of (x, v) during a period t is the backward trajectory over a period t of 
(x, —v). This raises the final proposition. 



Proposition 3.8. Let U be an open, bounded and C 1 domain in R . Then for all 
(x, v) inUxM. N we have existence and uniqueness of the characteristic (X t (x, v), V t (x, v 
and it is given by, for allt^O, 



X t (x,v) = X fin {t,X,-v) + (t-t fin (t,X,-v))Vf in {t,X, -v), 
V t (x,v) = -V fin (t,X,-v). 

Moreover, we have that V t (x,v) = Ot jXtV (v) with O t , x ,v an orthogonal transformation, 
and that for almost every (x, v) inlA x M. N we have the following 

(3.1) Vt ^ 0, (x, v) = (X t (X t (x, -v), -V t {x, -v)), V t (X t (x, -v), -V t (x, -v))). 



Proof of Proposition |3.8[ By construction we have that 

Ot,x,v — T^Xfin(t,x,v) ° • • • ° Tt xi (t,x,v)- 



It remains to show only the last equation (3.1 ), but this is direct this the backward 
trajectory of (x,v) is the forward trajectory of (x,—v). However, we can reach a 
point on D, stop after a time t\ and so the forward trajectory of that point during a 
time t > t\ does not come back to the original point (since we stayed in VL stop for a 
period t — t\). 

However, the set of points that reach VL stop belongs to the set of points that bounce 
infinitely many times in a finite time and this set is of measure zero (see Proposition 



3.4). □ 



3.3. Existence and uniqueness of solution to (2.5) — (2.7) 



3.4. Proof of uniqueness. The uniqueness of a solution with uq in C\ T\L 2 X v comes 
directly from the fact that we have a preserved quantity through time, thanks to 
the specular reflection property. Indeed, let us assume that u is a solution to our 
free transport equation satisfying specular boundary condition and the initial value 
problem uq. Then, a mere integration by part gives us 



vt ^o, iKt, v )ii; 



Nl 



which directly raises the uniqueness of a solution, since the transport equation (2.5) 
is linear. 
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3.4.1. Construction of the solution. It remains to construct a function u that will 
be constant on the characteristic trajectories and check that we indeed obtain a 
function that is differentiable in t and x which satisfies the transport equation. The 



first point of Remark 3.7 gives us the answer as we expect the following behaviour 

u(t, x, v) = u(t — ti(x, v),Xi(x, v),vi(x, v)) = ■ ■ ■ = u(t — tk(x, v),Xk(x, v),Vk(x, v )), 

up to the point where there is no more rebounds in a time t and then we go straight 
for the remaining period of time. 

Thus, we define: V(t,x,v) G R+ x U x R N , 

U(t, X, V) = U (X fin (t, X, V) - (t - t fin (t, X, V))v fin (t, X, v),Vf in (t, X, v)) . 

3.4.2. Boundary and initial conditions. First of all, u satisfies the initial condition 



(2.6) as n(0,x, v) = (since t min (x } v) ^ 0). 
u also satifies the specular boundary condition ( |2.7 ). Indeed, if (x,v) is in A + 



then either v ■ n{x) = and the result is obvious since lZ x (v) — v, or v ■ n{x) > 
and thus (x,lZ x (v)) belongs to Vt rebounds so t min (x,lZ x (v)) = (Proposition 3.3[ ). An 
easy induction raises 

x k (x,v) = x k +i{x,TZ x (v)), v k (x,v) =v k+ i(x,K x (v)), t k (x,v) = t k +i(x, TZ x (v )), 

for all k in N. 

The last equality gives us that n(t, x, v) = n(t, x, lZ x {v)) — 1 and therefore, combined 
with the two other equalities, 

X fin (t,X,v) = Xf in (t,X,TZ x (v)), V fin (t,X,v) = V fin (t, X, U x {v)), 
tf in (t,X,v) = t fin (t,X,Tl x (v)), 

which leads to the specular reflection boundary condition. 

3.4.3. Time differentiability. Here we prove that u is differentiable in time on IR + . 
Let us fix (x,v) in U x R N . 

By construction, we know that n(t,x,v) is piecewise constant. Since (tk(x, u))& e N is 
strictly increasing up to the step where it takes the value +oo, for t k (x,v) < t < 
tk+i(x,v) we have that for all s E R such that t k (x,v) <t + s < tk+i(x,v), 

X fin (t,X,v) = Xf in (t + S,X,V), Vf in (t,X,v) = V fin (t + S,X,v), 

t fin (t,x,v) = t fin (t + s,x,v). 

Therefore, we have that 

u(t + s, x, v) — u(t, x, v) 
s 

_ UpjXfin ~(t + S- t fin )v fin ,V fin ) - U (x fin - (t - t fin )v fin ,V fin ) 

s 

-> -Vfin ■ (V X M ) {Xfin - it- tfin)Vfin,Vf in ), 

because Uq is C 1 in x. So u is differentiable at t if t in strictly between two times 
t k {x, v). So we have that u is differentiable at t and that its derivative is continuous 
(since Xfi n , Vfi n and tf in are continuous when x and v are fixed). 

In the case t = tk(x,v) we can use what we just proved to show that we have 
the existence of right (except for t = 0) and left limits of d t u(t,x,v) as t tends to 
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tk(x,v). Using the specular reflection boudary condition of uq and the fact that it 
is C 1 in x gives us the equality of the two limits. 

3.4.4. Space differentiability and solvability of the transport equation. Here we prove 
that u is differentiable in x in U, which follows directly from the time differentiability. 
Let us fix t in IR + and v in M , we shall study the differentiability of u(t, -,v) in the 
direction of v. 
U is open and so 

Vx e U, 35 > 0, Vs e [-5, 5], x + sv e U. 

Thanks to the inductive construction, one find easily that 

u(t, x + sv, v) = u(t — s, x, v). 

Therefore, since u is time differentiable, we have that u(t,-,v) admits a directional 
derivative in the direction of v and that 

D x (v)(u)(t,x,v) = -d t u(t,x,v). 



4. MAXWELLIAN LOWER BOUND IN THE CUTOFF CASE: PROOF OF THEOREM 

In this section we are going to prove a Maxwellian lower bound for a solution to 



the Boltzmann equation (1.1) in the case where the collision kernel satisfies a cutoff 
property. 

The strategy to tackle this result follows the main idea used in [10] and [12] which 
relies on finding an "upheaval point" (a first minoration uniform in time and space 
but localised in velocity) and spreading this bound, thanks to the spreading property 
of the Q + operator, in order to include for larger and larger velocities. 

We gather here two lemmas, proven in [TU], that we will use a lot in this section. 
We remind the reader that we are using Grad's splitting (2.8). Let us first give an 
L°° bound on the loss term (Corollary 2.2 in [TP]). 



Lemma 4.1. Let g be a measurable function on 



Then 



Vv e 



\L[g)(v)\<Ci(v) 



where Cg is defined by: 



1) If <& satisfies (1.4) with 7^0 or if $ satisfies (1.5), then 



= est nfeC$e q . 



(2) //$ satisfies ([L4} with 7 e (-iV,0), then 

C L g = est n b C* [e g + $ , p > N/(N + 7). 

Then, our strategy is to spread the properties of the "upheaval point" thanks to 
a spreading property of Q + which is given by the following lemma (Lemma 2.4 in 
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Lemma 4.2. Let B = $6 be a collision kernel satisfying (1.3), with $ satisfying 
(1.4) or (1.5) and 6 satisfying (1.6) z/ ^ 0. T/ien for any v G K^, < r ^ R, 



£ G (0, 1), we have 

Q + {lB(v,R),lB(v,r)) ^ CSt l b C<i>r N ~ 3 R 3+1 Z^' 1 l B f v> ^/ r 2 +R 2 ( 1 ■ 

As a consequence in the particular quadratic case 5 = r = R, we obtain 



Q + (1 B (v,8),1b{v,6)) > est l b c^5 N+J ^ 2 1 \, 



5^(1-0)' 



/or any t> G R and £ G (0, 1). 



The case of the torus, studied in pU], indicates us that without rebound the 
expected minoration is created after time t = as quickly as one wants. Therefore 
we expect the same kind of bound to arise on each characteristic trajectory before 
its first rebound. However, in the case of a bounded domain, rebounds against the 
boundary can occur very close to the time t — and a rebound preserves only the 
norm of the velocity. Therefore, we will fail finding a uniformly (in space) small time 
where a uniform bound arises. Nevertheless, the convexity and the smoothness of 
the domain implies that grazing collisions against the boundary do not change the 
velocity very much. 

Thus our strategy will be split in three parts, which are the next three sections. 
The first step will be to partition the position and velocity spaces so that we have 
an immediate appearance of an "upheaval point" in each of those partitions. The 
second one is to obtain a uniform lower bound which will depend only on the norm of 
the velocity. Then the final part will use the standard spreading method used in [TO] 
and [12J which will allow us to deal with large velocities and derive the exponential 
lower bound uniformly. 

4.1. Partition of the phase space and first localised lower bounds. In this 



section we use the continuity of / together with the conservation laws (2.4) to 
obtain a point in the phase space where / is strictly positive. Then, thanks to 



the continuity of /, its Duhamel representation (2.9) and the spreading property of 
the Q + operator (Lemma 4.2) we extend this positivity to high velocities at that 
particular point (Lemma 4.3). Finally, the free transport part of the solution / will 
raise the immediate appearance of the localised lower bounds (Proposition 4.4). 

Moreover we define constants that we will use in the next two subsections in order 
to have a uniform lower bound. 

The study of the characteristic trajectories has already been done in Section [3] so 
we just define here new shorthand notations. 

For x in U, v in IR d and s, t ^ we denote the point at time s of the forward 
characteristic passing through (x, v) at time t by 

X S:t {x, v) = X s (X t (x, -v), -V t (x, -v)) 
V S!t (x, v) = V s (X t (x, -v), -V t (x, -«)), 



which has been derived from (3.1). 



We start by the strict positivity of our function at one point for all velocities: 
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Lemma 4.3. Let f be the mild solution of the Boltzmann equation described in 
Theorem 12.31 

Then there exists {x\, V\) inU xM, d and A > such that for all nGN, there exists 
r n > such that for all t in [0, A], there exists a n (t) > such that 

Vr G 5 (x h , Vf G R d , f(t,x,v) ^ a n (t)l 

with ao > independent of t and the induction formulae 



a n+1 (t) = C, 



N+<y / -min(£,A/(2"+ 2 r n ) 



4.1 



and 



where Cq = cstl^c® is defined in 4.2 and Cl = cstn^C^Ef (or Cl = cstnbCq,(Ef + 
L P j)) is defined in 



A 3v/ 2 
r = A, r n+1 = — — r r> 



Proof of Lemma |4T3j The proof is an induction on n. 

Step 1: Initialization. We recall the conservation laws satisfied by a solution 
to the Boltzmann equation, (2.4), 

V* G K + , 




f(t, x, v) dxdv = M, 




\v\ f(t, x, v) dxdv = E. 



Since U is bounded, and so is included in, say, 5(0, Rx), we also have that 



Vt G R 




(\x\ 2 + \v\ 2 ) f(t,x,v) dxdv ^a = MR 2 X + E < +oo. 



Therefore if we take t — and R m in = -\/2a/M, we have the following 



/ / 



M 



fo(x,v)dv ^ — > 

) J B(0,R min ) z 

Therefore we have that there exists xi in U and t>i in 5(0, R m in) such that 

r , N M 

4Vol(5(0,5 min )) i 

The first step of the induction is then due to the continuity of / at (0, Xx,Vi). 
Indeed, there exists A > such that 



Vt G [0, A], Vr G B(xi, A), Vv G B(vi, A), f(t,x,v) ^ 



M 



8Vol(5(0, i^in)) 2 



and we define A = min(5 r , <5x, £ 



Step 2: Proof of the induction. We consider the case where the Proposition 
is true at n. 

We take x in 5(0, A/2 n+1 ) and v in 5(0, + 2r n ) and t in [0, A]. 

We use the fact that / is a mild solution to write f(t,x,v) under its Duhamel 
form (2.9). The control we have on the L operator, Lemma 4.1| allows us to bound 
from above the integral term (the first term is positive). Moreover, this bound on L 
is independent on t, x and v since it only depends on an upper bound on the energy 
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e f(t,x,-) (and its local LP norm l^u x .\) which is uniformly bounded by Ef (and by 
Iff). This raises, for r n (t) = mm{t, A/(2 n+2 r„)) 



(4.1) 

f(t,x,v) ^ 







where Cl = est n b C^Ef (or Cl = est n b Cq>{Ef + see Lemma 4.1, and we used 
\\V Sit (x,v)\\ = \\v\\^2r n . 



Besides, we have that B(xi, A) cW and also 
A 



Vs G 



0. 



2 n+2 r 



Vw* G -8(0, || wi|| + 2r n ), ||xi — (x + su*)|| < 



A 



which, by definition of the characteristics (see Section 3.2), yields 



Vs G 



0. 



A 



2 n+2 r 



Wv, G B(0,||v a || + 2r„ 



= x + sv* G -B (x u 



Therefore, by calling the integration parametre in the operator Q + we can 



apply the induction property to f(s,X St t(x,v),v*) which raises, in (4.1) 



f(t,x,v) ^ 



sC L (\\v 1 \\+2r n )'r r n 2 



Applying the spreading property of Q + , Lemma 4.2, with £ = 1/4 gives us the 
expected result for the step n + 1 since B(vi,r n+ i) C B(0, \\vi\\ + 2r n ). □ 



Then we have all the tools to prove the next proposition which is the immediate 
appearance of localised "upheaval points". 



Proposition 4.4. Let f be the mild solution of the Boltzmann equation described 
in Theorem \2.3\ 

Then there exists A > such that for all < To ^ A, there exists St{to), o~x{t~o), 
3v{t~o), Rmin{T~o), ao( r o) > such that for all for all N in N there exists Nx in N* 
and x\, . . . , %n x inU and v%, . . . , v n x ^ n B(0, R m in{jo)) and 

• Uc (J B(x i ,8 x (r )/2 N ); 

• Vi G [r ,*r(r )],Vz G B(x h S x (r )), G JR d , f(t,x,v) ^ a (^o)lBK,<5 v (r ))(f ) 

5 (fi,5v(r )) C 5(0, -R mi „(r )). 

Proof of Proposition |4.4[ We are going to use the free transport part of the Duhamel 



form of / (2.9), to create localised lower bounds out of Lemma 4.3 



We take < To ^ A, where A is defined in Lemma 4J3 
U is bounded so let us define 

d u = sup {\\x — x'\\ , x, x' G U} 

and take n big enough such that ||t>i|| +r n ^ 2du/ro and define -R m m(To) = \\ v i\\ +r n - 
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(4.2) 



Thanks to Lemma 4.3| apply to this particular n we have that 

To 



Vt G 



2 ' 



A 



, Vx G B(x u A/2"), f(t,x,v) ^ a n (^) l B( „ l>rn 
where we noticed that a n (t) is an increasing function. 
Now we define 



m 

- 1 e 2 



2^ (I) 

We remark that / is continuous on the compact [r , A] x U x £>(0, -R m j n (ro)) and 
hence uniformly continuous. Therefore it exists S' t (to), S' x (to), <V( r o) > such that 

V|t - t'\ < ^(r ), V ||x - x'K *x fa), V - < <^(r ), 



(4.3) 



|/(t,x,v) - f(lf,a/,i/)\ ^ a (r ). 



We conclude our definition by taking 

5 t (t ) = min(A,r + S' t (t )), 
5 x (t ) = mm(5' x (r ),A/2 n ), 
S v (t ) = min(5y(r ),r n ). 

Finally, we take N G N and notice that U is compact so there exists x%, . . . ,xj^ x 
in U such that U C |J 5 (x^, 5x(to)/2 7V ). Moreover, we construct them such 



that Xi is the one defined in Lemma |4.3[ 

We then take V\ to be the one defined in Lemma |4.3| and we define 

2 

Vz G {2, . . . ,N X \, v i = —(x i -x 1 ) 



such that we directly have, because U is convex, 



(4.4) 



X TQ /2,Tx ) (Xi, Vi) — Xi, 
V T0 /2,T (Xi,Vi) = Vi, 



Now all the constants have been defined we can conclude by taking i in {2, ... , Nx} 
and using that / is a mild solution of the Boltzmann equation and therefore we can 
write it under the Duhamel form (2.9). As in the proof of Lemma 4.3 we can 



control the L operator appearing in the first term in the right-hand side of (2.9) 



(corresponding to the free transport). Thus, we use the duhamel form (2.9) between 



To and tq/2 and we combine it with (4.4). This yields 



f(r ,Xi,Vi) ^ f[-^,x u vi)e - 



T ^ -^C^ixt-xx))^ 



^ 2a (T )l B 



{vi,rn)\ V i)i 



where we used (4.2) for the second inequality. We recall here that Vi belongs to 



B(0,R m i n )(T )) and that 5(0, i? min (r )) C B{v x ,r n ) and therefore 
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(4.5) f(ro,Xi,Vi) ^ 2a (r ). 



Finally, combining (4.5) with the uniform continuity of /, (4.3) we have that for 



all t in [tq, 5t(tq)], x in B(xi, 5x(tq)) an d v in B(v i} 5y(r )), 

f(t,x,v) ^ a (7"o)- 



□ 



Remark 4.5. 77ms Zas£ Proposition tells us that localised lower bounds appear im- 
mediately, that is to say after all time tq > 0. The exponential lower bound we 
expect will appear immediately after those initial localised lower bounds, that is to 
say for all t\ > t . Therefore, to shorten notation and lighten our presentation, 
we are going to study the case of solution to the Boltzmann equation which satisfies 



Proposition 4.4 up to r = 0. Then we will immediatly create the exponential lower 



bound after and apply this result to F(t, x, v) = f(t + r , x, v). 

4.2. A lower bound only depending on the norm of the velocity: proof 
of Proposition |2.4[ As noticed in the introduction, the spreading property of the 
bilinear operator Q + cannot be used (at least uniformly in time and space) when 
we are really close to the boundary due to the lack of control over the rebounds. 
However, if we have a lower bound only depending on the norm of the velocity then 
the latter bound will not take into account rebounds as they preserve the norm. 
Therefore we will be able to spread this minoration up to an exponential one. 

This section is dedicated to the creation of such a uniform lower bound depending 
only on the norm of the velocity. In order to do so we restrain the problem without 
taking into account large velocities and divide the study in the two cases whether 
the trajectory stays close to the boundary or not. In both cases we will start from 



the localised "upheaval points" constructed in Section [4TT] and spread them to the 
point where one gets a lower bound depending only on the norm of the velocity. 

The next subsections tackle each of these points. We first study the case when 
a characteristic reaches a point far from the boundary and finally we focus on the 
case of grazing characteristics. We recall here that in the next two subsection we 



consider functions that satisfies the localised lower bounds, Proposition 4.4, at the 



time t = 0. Thus, we fix St, 5x, Sy, Rmin and ao to be the ones described in 



Proposition 4.4 at time r = 0. 



The result we will derive out of those studies is Proposition 2.4 and for now on, 
dependencies on physical observables of / (Ef and If^) will be mentionned when 
constructing our constants and we will not explicitely be written everytime. 

4.2.1. Characteristics passing by a point far from the boundary. In this section we 



manage to spread the lower bounds created in Proposition 4.4 (with r = 0) up to a 
ball in velocity centred at zero as long as the trajectory we look at reaches a point 
far enough from the boundary. 

First, we pick N in N* and we can cover U with [J 1<i<Nx B(xi, 5 X /2 N ) thanks to 
Proposition |4.4 



Then for Z ^ we define 
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(4.6) Ui = {xeU: d(x, dU) ^ 1} , 

where d(x, dU) is the distance from x to the boundary of hi. 

We now take R > and we define two sequences in M + by induction, for all r ^ 
and I ^ 0, 



(4.7) 



(4- 



and 




Cr, r " ' ' L-p-rCiW n 2 ( l - t) 



where Cq and Cl have been defined in Lemma 4.3 



We express the spreading of the lower bound by the following proposition. 



Proposition 4.6. Let f be the mild solution of the Boltzmann equation described 



in Theorem \2.3\ and suppose that f satisfies Proposition ^. 4| with ro = 0. 

Consider < r ^ 5t and N in N. We consider (a^ien,...,.^} an d ( v i)ie{i,...,N x } 

described in Proposition \4A\ with To = 0. 

Then for all n in {0, . . . , iV} we have that the following holds: for all < I ^ 5x, 
for all R > such that l/R < r, for all t in [I / (2 n R) , t] , and for all x G U and 
v G B(0,R), if there exists t\ G [0,t — l/(2 n R)} such that X tu t(x,v) belongs to 
UinB(xi,8 x /2 n ) then 

f(t,x,v) ^ a n (l,T)l B(v ^ rn) (V tut (x,v)), 



where (r n ) and (a n ) are defined by (4.7)-(4.8) 



Proof of Proposition 4.6 This Proposition will be proved by induction on n. 



Step 1: Initialization. The initialisation is simply Proposition 4.4 and the first 



term in the Duhamel formula (2.9) starting at r 



Indeed, we use the definition of / being a mild solution to write f(t,x,v) under 
its Duhamel form (2.9) starting at ti where both parts are positive. The control we 



have on the L operator, Lemma 4.1, allows us to bound from above the first term. 



Moreover, this bound on L is independent on x and v (see proof of Lemma 4.3). 
This raises 



(4.9) 



f{t, x, v) > e-e-'O^nr f ( tl ,X tl>t (x, v), V tl>t (x, v)) . 



Finally, Proposition 4.4 applied to f(ti, X tl>t (x, v), V^ li4 (x, v )) gives us the property 
for n = 0. 

Step 2: Proof of the induction. We consider the case where the Proposition 
is true at n. 

We take I G (0, S x ], t G [l/(2 n+1 R), r], x G U and v G 5(0, R). 

We suppose now that there exists t\ G [0, t — l/(2 n+1 R)] such that X tl t{x,v) G 
lA^B^Sx/Z^ 1 ). 
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Following what we did at the first step of the induction but on the integral part 



of the Duhamel formula (2.9) we reach 



(4.10) f(t,x,v)> 

e -C L r(RV + ( f» ~ -Q+[ f ^ Xs!t {x,v),-)J(s,X S!t {x,v),-)} ds I (V tl>t (x,v)) 



tl+ 2 „ + 2 R 



The goal is now to apply the induction at rank n to the triplet (s, X s j(x, v), v*), 
where t>* is the integration parametre inside the Q + operator, with ^ R. 

One easily shows that X Sit (x,v) = X tl j(x,v) + (s — ti)V tlt t(x,v), for s in [t\,ti + 
l/(2R)], and therefore we have that 

(4.11) \\X h>t (x,v)-X Sit (x,v)\\^ ' 



2n+2 ' 

and so that X s j(x, v) belongs to 

Finally, we have to find a point on the characteristic trajectory of (s, X s j(x, v),v*) 
that is in Uy for some This is achieved at the time t\ (see FigJ4l) . 




FIGURE 4. Study of (s,X Stt {x, v),v*) far from the boundary 
Indeed, we have s in ft + l/(2 n+3 R), t x + l/{2 n+2 R)\ so, for < R 



(4.12) 



W e [ti,s], \\X s j(x,v) - (X S}t (x,v) - (s-s')v*)\\ ^ 



2 n + 2 ' 



This gives us the characteristics trajectory backward starting from s, since iX s j(x, v) — 
[s — s')v* remains in W, and therefore 



Vs'eft,*] 



X 8 ',s {X s , t {x, v), u*) = X s>t (x, v) - (s - s')u* 
V s >, s (X S:t (x,v),v*) =v*. 
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To conclude we just need to gather the upper bounds we found about the trajecto- 
ries reaching (X Sjt (x, v), v*) in a time s in [ti + l/(2 n+3 R),ti + 1 / (2 n+2 R)] ) equations 
CTIj) and p7i2j)' 



\X tlt (x,v) - X tl>a (X at (x,v),v*)\\ < 



2 n+l 



We have that X tl j(x,v) belongs to Ui PI B(xi, S x /(2 n+1 )) and therefore we have 
that for all s in [tj +'//(2 n+3 J R), *i +l/{2 n+2 R)}, X tl s (X s t (x, v),v m ) belongs to U l/2 n 
B(xi,6 x /F). 

Finally, if s belongs to [t x + //(2 n+3 J R), ^ + l/(2 n+2 R)] we have that (1/8) /(2 n R) < 
s ^ r and t x is in [0, s - (l/8)/(2 n R)]. 

We can therefore apply the induction assumption for I' = 1/8 inside the Q + 
operator in ( 4.10[ ), recalling that Vt ltS (X Stt (x, v), v*) = v*. 



f(t,x,v) ^ a r , 



r e 



-C l t{R)^ 



*l+ 2 n + 2 H 




<2 + [lB(« 4 ,r„),lB(t» i ,r n )] ^ (Vfc.tfo «)) 



2 n + :i fi 



Applying the spreading property of Q + , Lemma 4.2, with £ = 1/4 gives us the 
expected result for the step n + 1. □ 

One easily notices that (r n ) ne pj is a strictly increasing sequence. Moreover, for all 
N in N we have that for all 1 ^ % ^ N x , Vi belongs to B(0,R m i n ). Therefore, by 
taking N big enough (bigger than N\ say) we have that 

Vi G {1, . . . , N x }, B(0, 2R mm ) C B(v h r N ). 
This remark leads directly to the following corollary which stands for Proposition 



2.4 in the case when a point on the trajectory is far from the boundary of 14. 



Corollary 4.7. Let f be the mild solution of the Boltzmann equation described in 



Theorem 2.3 and suppose that f satisfies Proposition 4.4 with Tq = 



Let At be in (0,5t) and take T\ in (0, At]- 

Then for all < I ^ 5 x , there exists o(1,t\,At) > and < a(l,Ti, At) < T\ 
such that for all t in [t\,At], and every (x,v) in U x M. d : if there exists t\ G 
[0, t — a (I, Ti, At)] such that X tl j(x, v) belongs to Ui then 

f(t,x,v) ^ a(l,Ti,A T )l B (0,2R min )(v). 



Proop of Corollary 4.7 This is a direct consequence of Proposition 4.6 



Indeed, take < / < 5 X , < n < A T and R = R(A T ) > such that R ^ 3R min 
and l/R < A T . Then take N 2 ^ N x big enough such that l/(2 N2 R) < t x . We 
emphasize here that N 2 is defined thanks to t% so we write N 2 (ti). 

Now apply Proposition 4.6 with iV = N 2 (t{) and for t in [r^A^]. We obtain 



exactly Corollary 4.7 (since B(0,2R min ) C B(v i7 r x ) for all z and i? ^ 3R min ) with 



a(/,ri, A T ) = a A r 2 ( Tl) (/, A T ) and a(/,Ti,A T ) 



2^f2(A r )' 
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and the fact that |J B (xi,5 x /2 N ) covers U. □ 

4.2.2. Grazing trajectories. We now turn to the case when the characterictic tra- 
jectory never escapes a small distance from the boundary of our convex domain 
U. 

Intuitively, by considering the case where U is a circle, one can see that such a 
behaviour is possible only either the angles of collisions with the boundary remain 
really small (which corresponds in higher dimension to the scalar product of the 
velocity with the outside normal being close to zero) or the angle is important 
but the norm of the velocity or the time of motion is small. Thus, by using the 
spreading property of the Q + operator we may be able to create bigger and bigger 
balls in between two rebounds against the boundary because the latters should not 
change the velocity too much. 

The study of grazing collisions will follow this intuition. First of all we will prove a 
geometric lemma dealing with the fact that if the velocities are bounded from below 
and above, then in small times, the possibility for a trajectory to stay very close to 
the boundary implies that the velocity do not change a lot along time. Then we will 
spread a lower bound, in the same spirit as in the last subsection, up to the point 
when this lower bound covers a centred ball in velocity. Notice that the geometric 
propery imposes us to work with velocities of which the norm is bounded from below 
and so we shall have to take into account the speed of the spreading. 

Geometric study of grazing trajectories. 

The key point of the study of grazing collisions is the following geometric lemma. 
We emphasize here that this is the only part of the article where we need the fact 
that U is convex and C 2 . 

Proposition 4.8. Let U be an open convex bounded C 2 domain in M. d and take 
< v m < v M - 

For all e > there exists t £ {vM) such that for all < r 2 ^ t e (vM) there exists 
l £ (t) mi T 2 ) > such that for all x inU and all v in M with v m ^ \\v\\ ^ Vm, 

(Vs G [0,r 2 ], X s {x,v) iUi e{Vm , T2) ) => (Vs G [0,t e (v M )], \\V s {x, v) - v\\ < e) . 
Furthermore, l e (v m ,-) is an increasing function. 

The proof of this proposition is quite technical and completely independent of the 
study of a lower bound for solutions to the Boltzmann equation. Therefore, we leave 
it to Appendix [A] for the sake of clearness. 

Now we understand how grazing trajectories behave geometrically we can turn to 
their effects onto the spreading property of the Boltzmann Q + operator. 

Spreading effect along grazing trajectories. 

In order to use the geometrical behaviour of grazing characteristic trajectories, the 
main issue is to consider velocities bounded from below. However, in the same time, 
we would like to spread a lower bound up to ball centred at which seems to fail 
the possibility of having a lower bound on the norm of velocities. We shall overcome 
this problem thanks to the flexibility of the spreading property of the Q + operator, 
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Lemma 4.2 which allows us to extend the radius of the ball from up to a/2 times 
the initial radius. 

The idea is to spread the initial lower bound by induction as long as the origin 
is strictly outside. This will allow us to use the geometrical property of grazing 
characteristics. Finally, a last iteration of the spreading property, not requiring any 
a priori knowledge on characteristics, will include in the lower bound. 

In Corollary 4.7 we can fix a special time T\ of crossing the frontier of some Ui 



and this allows us to derive a lower bound for our function in this special case. The 



second case of grazing trajectories is dealt with Proposition |4.8| where we can find 
an I for U\ to control the evolution of the velocity. Our goal now will be to find all 
the constants that are still free and to finally find a time of collision small enough 
that it will remain the same during all the iteration scheme. 



We now fix all the constants that remain to be fixed in Corollary 4.7 thanks to 



Proposition 4.8 



We first build 



(4.13) A r = min (S T , t Sv/i (3R min )) . 

Second of all we define, for f in (0, 1/2 - 3/(8^2)}, 



(4.14) 



r. 



ro(0 =5 V 

l+1 (0 =V2(1-£K(0 



4 ' 



We have that (r n (l/2 - 3/(8^2))) is a strictly increasing sequence. Therefore, 
it exists N max such that 



Now we fix N in N* greater than N max . With this iV and Proposition |4.4| at 
t = 0, we construct v%, . . . , vn x - 

For i in {1, ... , N} we take £® in (0, 1/4 -3/(16^2)] and we define N max (i) to be 
such that 0£B («i,r n (f«)) for all n < N max (i) and G B (v h r Nmax{ ^)). And 
we can in fact take £w SU ch that G Int (_B (t> i? rN max (i)(^))) ■ 

Therefore we have that for alH in {1, . . . , Nx}, 

* = INI-^ m- . ( 4)-i(^ i) )>o l 

which is strictly positive if and only if N max (i) > 0. We consider 



(4.15) v m = min 8 { > 0}. 

i£{l,...,N x } 



We can now define: 
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(4.16) 



R(t) = max I 3R 



25 



x 



(4.17) 
(4.18) 

Finally, we define l(r) 



n(r) = r - 



25 



x 



(4.19) 
(4.20) 



VO < r ^ A T , 



a t 



T2 T 



fl(r)' 
«(/(r),n(r),AT). 



J2(r) 



Z(r) = min (<J X , Zj y/4 (u m , r 2 (r))) . 
We also build up the following sequence, where R, I and T\ depend on r, 
(4.21) 



age 



-(A t ~t)C l {R)^ 



b^r, A T ) = min (^^(^^^-^e-^^fo^^ A T ) 2 ; a 2 (Z, r u A T ) 
£W was defined above and a(l, r, At) was defined in Corollary 



4.7 



We have the tools to state the next Proposition which is the complement of Propo- 



sition 4.6 in the case when the trajectory stays close to the boundary. We remind 



the reader that < a(r) < tAt). 



Proposition 4.9. Let f be the mild solution of the Boltzmann equation described 



in Theorem 2^ and suppose that f satisfies Proposition AA with tq = 0. 
Consider < r ^ At and take i in {1, . . . , Nx} such that N max (i) > 1. 
For all n in {0, ... , N max (i) — 1} we have that for all t in [r — 5x/{2 u R(t)), At], all 
x in B(xi,d~x/2 n ) and all v in B(0, R(t)) , if 

Vs e [0,t-a{r)), X Syt {x, v) W 1(t) 

then 

f(t,x,v) ^ b^(r,A T )l B ^ rnm ^(v), 



all the constants being defined in (4.13), (4.14), (4.20), (4.16), (4.17), (4.18) and (4.21). 



Proof of Proposition 4.9 We are going to use the same kind of induction we used to 
prove Proposition 4.6 So we start by fixing i such that N max {i) > 1. 



Step 1: Initialization. The initialisation is simply Proposition 4.4 and the first 



term in the Duhamel formula (2.9) starting at r, with the control from above on L 
thanks to Lemma [4.11 



Stef 2: Proof of the induction. We consider the case where the Proposition 
is true at n ^ N max {i) — 2. 

We take t in [r - 5 x /(2 n+1 R(r)), A T ), x in B(x u S x /2 n+1 ) and all v in B(0,R(r)). 

We suppose now that for all s G [0, t — a(r)] we have that X Sjt (x,v) does not 
belongs to Uu T \. 
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To shorten notation we will skip the dependence in r of the constant. 
We use the definition of / being a mild solution to write f(t,x,v) under its 
Duhamel form (2.9) where both parts are positive. As in the proof of Proposi- 



tion 4.6 we control, uniformly on t, x and v, the L operator from above. This 



raises 



" 2 n + 2 R 



' 2 n + 1 R 



(4.22) 

f(t,x,v) ^ e-° LT{R)1 ' 

where we used ||V^ )t (a;, v) 
holds true thanks to the definition of (4.16) 

Sx 



Q + [f(s, X Sjt (x, v), •), f(s, X Sjt (x, v), •)] (V Stt (x, v)) ds, 



\v\\ ^ R. We also emphasize here that this inequality 



The goal is now to apply the induction at rank n to the triplet (s,X S)t (x, v),v* 
where v* is the integration parameter inside the Q + operator, with ||t>*|| ^ R. 
We notice first that for all s in [t - S x /(2 n+1 R), t - 5 x /(2 n+2 R)} 



\xi-X S)t (x,v)\\ < 



5 



X 



2n+l 

Sx 



x - X Sit (x,v) 
Sx 



, At- s)R ^ 



so that for all s in [t-S x /{2 n+1 R), t-5 x /{2 n+2 R)}, X s , t (x, v) belongs to B{x h S x /2 n ). 
We also remark that 



_ S x S x 



2n+lR' 2 n + 2 R 



c 



5x A 

t , At 

2 n R' 



Now we have two different cases to consider for (X s r )S (X S) t(x, v), u*)) a / e[o,s-a]- 
Either for some s' in [0, s — a], X s r }S (X S} t(x, v), u*) belongs to Ui and then we can 



apply Corollary 4.7 



f(s,X Stt (x,v),v m ) ^ a(l,n, A T )l B (o i3Rmin )(v m ) 
(4-23) > b^(r,A T )l B{vurnm) (v), 

since V{ is in B(0, R m i n ). 

Or for all s' in [0, s — a] C [0,r 2 ], X s > jS (X Sjt (x,v),v*) does not belong to Ui and 
then we can apply our induction property at rank n and we reach the same lower 
bound pi. 



Plugging (4.23) into (4.22) raises, thanks to the spreading property of Q + , Lemma 
42] with £ = ■ fW, 
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(4.24) f{t,x,v)> 



' 2 n + 2 R 



To conclude we use the fact that for all s in [0, t — a] we have that X Sji (x, v) does 
not belong to IA\ and that t — a > r%. Moreover, n + 1 ^ N max (i) — 1 and so if v 
belongs to B (i>i,r n (£W)) we have that v m ^ ||u||. We apply Proposition |4.8[ raising 



Vs g 



?X 



\v-V Stt (x,v)\\ < 



5 



V 



2 n+l jR ' 2 n + 2 R_ 

Therefore, if v belongs to B (yi,r n+ i (fW)) we have that V s . i(x,v) belongs to 
B(vi, y/2(l - £ (i) )r„(£ (i) )) for all s in [t - 5 x /(2 n+1 R), t - 5 x /(2 n+2 R)]. 

Therefore if v belongs to B (v i, ?"n+i(£^)) we can compute explicitely (4.24) and 
obtain the expected induction. □ 



Thanks to Proposition |4.9[ we can build, for all x and all v, a lower bound that 
will contain in its interior after another use of the spreading property of the Q + 



operator. The next Corollary is the complement of Corollary 4.7 



Corollary 4.10. Let f be the mild solution of the Boltzmann equation described in 
Theorem \2.3 and suppose that f satisfies Proposition 4A with tq = 0. 
Let At be defined by (4.13). 



There exists ry > such that for all r G (0, At] there exists b(r) > such that for 
all t in [r, A T ] 



If for a(r) and l(r) being defined by ( |4.18[ ) - Q4.20 ), 

Vs G [0,t — a(r)], X s>t (x,v) £ U l{r) . 

Then 



f{t,X,v) ^ 6(r)lB(0,r(T)(«)- 



Proof of Corollary (4.10). As we said before stating our result, we are going to use 



the spreading property of Q + one more time. By definition of N max (i), for all , 



Vz G {1, . . .,N X }, G Int (B {v u r Nmax{t) (^))) 



We define 



r v = mm{r Nmax{i) (£ 



Will ; % 



e{l,...,N x }} 



which only depends on 5y and (fi)ie{i,...,7V x }- By construction we see that 



(4.25) 



Mi G {l,...,N x },B(0,r v ) C B {v t ,r Nmax{t) (^)) . 



Now we take r in (0, Ay] and we take t in [r, At], x in B(xi,8 x /2 Nx ) and v in 
B(0, R(t)) such that 
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Vs E [0,t-a{r)], X s>t (x, v) $ U 1{t) , 
We have that t is in [r - 6 x /{2 Nma " W _1 i2(r)), A T ] and x in B(x il 6 x /2 Nmm ®~ 1 )- 



By the same methods we reached (4.24), we raise for n = N max (i) 



(4.26) f(t,x,v)> 

C Q r^{^) N ' 2 -^ c ^\b®) 2 



t- 



2 n + 2 R 



1 B(* i ,^(l-{W)r„({(0)) (^(X.V)) ds. 



This time the conclusion is different because we cannot bound the velocity from 



below since our lower bound contains 0. However, (4.25) allows us to bound from 



below the integrand in (4.26) by a function depending only on the norm. Besides, 



IMI = ||^s,t(^, v)\\ along characteristic trajectories (see Proposition (3.8)). Thus we 
obtain the expected result by taking 



b(T)=wm\b<t> mamW] ie{l,...,N x }\. 



4.2.3. Proof of Proposition (2.4). By combining Corollary 4.7 and Corollary 4.10 



□ 



we 



can deal with any kind of characteristic trajectory and this is stated by the following 



lemma as long as our solution / satisfies Proposition 4.4 with r = 0, meaning that 
we have localised lower bound at time 0. 



Lemma 4.11. Let f be the mild solution of the Boltzmann equation described in 



Theorem 2.3 and suppose that f satisfies Proposition 4.4 with r = 0. 

There exists At > and ry > such that for all < r ^ Ay there exists a(r) and 



Vt E [t, A t ], a.e. (x,v) EU x 



f(t,x,v) ^ a(r)l B (o )Pv )(u). 



Proof of Lemma 4.1 1[ In Corollary 4.10 we constructed At and ry. 

We now take r in (0, At] and consider t in [r, At], (x, v) in U x R d where / is a 
mild solution of the Boltzmann equation. 



We remind the reader that l(r) and a(r) have been introduced in (4.20) and 



(4.18) 



Either (X Stt (x, v )) a e[o,t-a(r)] meets Ui{ T ) and then we use Corollary 4.7 
f(t,x,v) ^ a(Z(r),ri(r), A T )l B (o,rv)( v )- 



to get 



Or (X Sjt (x, u)) s e[o,t-a(T)] stays out of i^( T ) and then we use Corollary 4.10 

f(t,x,v) > &(t)1b(o,»v)( v )- 



to get 



We obtain Lemma 4.11 with a(r) = min (a(/(r), 7~i(t), At), &(t)). 



□ 
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The we have all the tools to prove Proposition 2A Indeed we can build a lower 
bound ony depending on the norm of the velocity immediately after having created 



the initial localised lower bounds in Proposition |4.4[ which appears immediately as 
well. 



Proof of Proposition 2.4 Let r be strictly positif and consider t such that t belongs 
to[r/2,r]. 



First case. First we suppose that / satisfies Proposition 4.4 with r = 0. 



We can compare t with A^ constructed in Lemma 4.11 



If t ^ At then we can apply the latter lemma and raise for almost every (x, v) in 
U x R d 



(4.27) 



f(t,x,v) >b[-\ l B (o, rv )(v) 



If t ^ Ax then we can use Duhamel formula ( |2.9 ) for / is a mild solution of the 
Boltzmann equation. We can then bound f(t,x,v) by its value at time Ay (as we 
did in the first step of the induction in the proof of Proposition 4.6 ) and use Lemma 



4.11 at At- This raises, for \\v\\ ^ ry, 



f(t,x,v) > f(AT,X AT , t (x,v),V AT , t ( Xj v))e-^ A ^ c ^ 
> b(A T )e-^- A ^ c ^ 1+ l B{0!rv) (V AT , t (x,v)) 
(4.28) ^ b{A T )e-^-^ c ^ 1+ l B{0trv) {v). 



We just have to take the minimum of the two lower bounds (4.27) and (4.28) to 



obtain Proposition 2.4 



Second case. We do not assume anymore that / satisfies Proposition 4^4 with 

T = 0. 



Thanks to Proposition 4.4 with tq = r/4 we have that 

Vt < 0, Vx G U, v G R d , F{t, x, v) = f(t + t , x, v) 
is a mild solution of the Boltzmann equation satisfying exactly the same bounds as 



/ in Theorem 2.3 and such that F has the property of Proposition 4.4 at (note 
that all the constants depend on r ). 

Hence, we can apply the first step for t' in [r/4, 3r/4] and F(t',x,v). This gives 
us the expected result for f(t, x, v) for t — t' + To in [r/2, r]. 

□ 
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4.3. Proof of Theorem 2.3[ As said in Section L2 the main difficulty was to 
create a lower bound depending only on the norm of the velocity. This has been 



achieved thanks to Proposition 2.4 The latter, consider as a starting point of an 
induction, can lead to exactly the same process developped by Mouhot in [Tu] , 



Section 3. Therefore we will just explain how to go from Proposition 2.4 to the 



complete Theorem 2.3, without detailing too much when unnecessary. 



First of all, by using the spreading property of the Q + operator once again we can 



grow the lower bound derived in Proposition 2.4 



Proposition 4.12. Let f be the mild solution of the Boltzmann equation described 
in Theorem \2.3\ 

For all t in (0, T), there exists Rq > such that 

T 



Vn e N, Vt G 



T 



2 n+1 

with the induction formulae 



a 2 , T \ T N+jtN/2+l 

a n +i(r) = est C e wV " and r n+1 = V2(l - £ n )r„, 



2n+l 

where (^ n )neN is any sequence in (0, 1) and r = ry, Oq{t) and C e only depend on 



r, E f (and Ly if $ satisfies (1.4) with 7 < 0) 



Indeed, we take the result in Proposition |2.4| to be the first step of our induction 
and then, for ninN and < r < T, the Duhamel form of /, with the usual controls 
we made above, raises 



f(t,x,v) ^ 



- 2 n + l 

fort in [r-r/2 n+2 ,r]. 

Using the induction hypothesis together with the spreading property of Q + (Lemma 



4.2 ) leads us, as in the proofs of Propositions 4.6 and 4.9 , to a bigger ball in velocity 



centred at 0. The only issue is to avoid the v -dependence in exp —C^it — s)(v ) 7+ 
which can easily be achieved as shown in the end of the proof of Proposition 3.2 in 

This is exactly the same result that pU], Proposition 3.2, but uniformly in x as we 
managed to derive Proposition 2.4 uniformly in space (which seems to be different 
with the result in [10] but is not as Mouhot implicitely assumed that the initial data 
was uniformly in space bounded from below). 

As in Lemma 3.3 in [TU], we can take an appropriate sequence (£ n )neN an d look 
at the asymptotic behaviour of (o n (r)) ngN . We obtain the following 



Vr > 0, 3p T , 9 T > 0, V(x, v) G U x 



f{t,x,v) > 



Pr 



_H1 
e 2« 
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Notice here again that the result is uniform in space, since the previous one was, 
and that the constants p T and 9 T only depend on r and the physical quantities asso 
dated to /. 



To conclude it remains to make the result uniform in time. As noticed in [10] . 
Lemma 3.5, the results we obtained so far do not depend on an explicit form of fo 
but just on uniform bounds and continuity that are satisfied at all times, positions 
and velocities. Therefore, we can do the same arguments but not starting at t — 
but at any time. So if we take r > and consider r ^ t < T we just have to make 



the proof starts at t — r to obtain Theorem 2.3 



5. Exponential lower bound in the non cutoff case: proof of 

Theorem 



In this section we prove the immediate appearance of an exponential lower bound 



for solutions to the Boltzmann equation (1.1) in the case of a collision kernel satis- 
fying the non cutoff property. 

The definition of being a mild solution in the case of a non cutoff collision kernel, 



Definition 2.5 and equation (2.11 ), shows that we are in fact dealing with an almost 
cutoff kernel to which we add a non locally integrable remainder. The strategy will 
mainly follow what we did in the case of a cutoff collisin kernel and to compare it, 
at each step, with the loss due to the remaining part. 

As in the last section, we shall first prove that solutions to the Boltzmann equation 
can be uniformly bounded from below by a lower bound only depending on the norm 
of the velocity and then use the proof given for the non cutoff case in [10] . Basically, 
we prove the immediate appearance of localised "upheaval points" and spread them 
up to the point where we reach a uniform lower bound that includes a ball in 
velocity centred at the origin. The spreading effect will be done both in the case 
where the trajectories reach a point far from the boundary and in the case of grazing 
trajectories. Then we will spread this lower bound on the norm of the velocity up 
to the exponential lower bound we expect. 

We gather here two lemmas, proved in JTU], which we shall use in this section. 
They control the L°°-norm of the linear operator S e and of the bilinear operator Ql. 



We first give a property satisfied by the linear operator S, (2.11 ), which is Corollary 
2.2 in [H]. 

Lemma 5.1. Let g be a measurable function on R . Then 

\JveR N , |%](«)|<C»^, 

where is defined by: 



(1) If satisfies (1.4) with 7^0 or if Q satisfies (1.5), then 



C g = est mbC<z>e g . 
(2) //$ satisfies (O) with 7 G (-iV,0), then 



C s g = est m b C*[e g + lP], p>N/(N + j). 

As announced above, we compare the lower bound created by the cutoff part of 
our kernel to the remaining part Q\. To do so we need to control its L°°-norm. This 
is achieved thanks to Lemma 2.5 in [101 . which we recall here. 
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Lemma 5.2. Let B = $6 be a collision kernel satisfying (1.3), with $ satisfying 



(1.4) or (1.5) and b satisfying (1.6) z/ G [0,2). Lei /, g be measurable functions 

on R d . 

Then 



[1) If $ satisfies (1.4) 2 + 7^0 or i/ $ satisfies (1.5), i/ien 
^ < est m 6 C ||g|| Li 



Vv e 



(2) J/ $ satisfies ( |L4| ) 2 + 7 < 0, t/ien 
Vt;GM d , < est m 6 C* 

withp> N/(N + 7 + 2). 



The strategy used for the cutoff case is reproduced here. That is to say, our first 
aim is to construct a first lower bound only depending on the norm of the velocity 
and this is achieved by studying separately characteristic trajectories going far from 
the boundary and the others staying close to it. 

5.1. A lower bound only depending on the norm of the velocity. In this 



section we prove the following proposition, which is exactly Proposition 2A but 
applied to non-cutoff collision kernels. 



Proposition 5.3. Let f be the mild solution of the Boltzmann equation described 
in Theorem \2.Q[ 

For all < r < T there exists ao( r ) > such that 

Vt G [r/2,r], V(x, v) G U x R d , f(t,x,v) ^ a (r)l B{0irv) (v), 



ry and ao(r) only depending on Ef, E'r, Wf (and If? i/$ satisfies (1.4) with"/ < 0) 



Proof of Proposition 5.3 As before, we would like to create localised "upheaval 



points" (as the ones created in Proposition 4.4) and then extend them. Both steps 
are done, as in the cutoff case, by induction along the characteristics. 
We have the following inequality 

(5.1) Qtif, f) + Ql(f, f) > Qtif, f) - \Ql(f, f) I , 



the definition of being a mild solution in the non-cutoff case, see Definition 2.5 , for 
all < e < Eq, 



(5.2) 

f(t, X t (x, v), V t (x, v)) = f (x, v)exp 



- / (L e + S E )[f(s,X s {x,v),-))(V s (x,v)) ds 



+ / exp 
'0 



(L £ + S £ ) [f(s , ,X s ,(x,v),-)](V sl (x,v))ds' 



(Qt + Ql) [f(s, X s (x, v),-), f(s, X s (x, v), -)}(V s (x, v)) ds. 
and the following bounds thanks to Lemmas 4T , |5.1 and |5.2| 



(5.3) 
and 



Lelf] < C f n b co(vy + , S £ [f] < C f m b Nco(vY 
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(5.4) 



\Ql(f,f)\^C f m b Nco(v)( 2+ ^ 



with Cf > is a constant depending on Ey, E'j, Wf (and L^ 7 if $ satisfies (1.4) 
with 7 < 0). 



The proof of Proposition 5.3 is divided into three different inductions dealt in the 



same way as in the proof of Proposition 2.4 One to create localised initial lower 



bounds (Lemma 4.3), another one to deal with non-grazing trajectories (Proposition 



4.6) and the final one for grazing trajectories (Proposition 4.9). Therefore, we will 



just point out below the only changes we need to make to make those inductions 
work in the non-cutoff case. 

In all the inductions in the cutoff case, the key point of the induction is to control 
at each step quantities of the form 



f(t,x,v) ^ 



4V 



(i) 



exp 



(L £ + S e ) [f(s\X s ,(x,v), -)}(V s ,(x,v)) ds' 



(Qt + Ql) [/(*, X a {x, v), •), f(s, X s (x, v), -)}(V s (x, v)) ds, 

(tn ) n es, (in')neN being defined differently for grazing and non-grazing trajectories 
(see proofs of Propositions 4.6 and 4.9). 



Exactly like in those previous induction, and using (5.1), (5.2) and (5.3) — (5.4) 
if f{t,x,v) ^ a n l B (v,r n ) then 



„4 2 ) 

f(t,x,v) > e~ c f^ (a 2 n Qt[l B (v,r n ),lB(v,r n )} -C f m b Nco(R)^ + ) (V s (x,v))ds, 
J tit' v 7 



which leads to 



(5.5) f{t,x,v)> 



,(2) 



e -Cf(R) 



(i) 



a^cst hcoc^r^Zn 2 ^(^^2(1-60) ~ C f m b Noo(R)^ + ^ j (V s (x,v)) ds 



(2+7)^ 



by using the spreading property of Qf (see Lemma 4.2) and using the shorthand 
notation Cj{R) = Cf(n b co + m b Nco){R)' y . 

To conclude we just notice that 

l b co l b 

and that 

(5.6) - b ° " 6n 

if v belongs to (0, 2) and 
(5.7) 

for v = 0. 



-e 2 ~ u 



n b co ~ — £ , rrihNco ~ 



b 



n b co ~ b | logs | , m b Nco - 

e e-s>0 e e->0 I — V 
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Thus, at each step n of the inductions we just have to redo the proofs done in the 
cutoff case and choose e = e n small enough to be such that 

(5.8) C f m b Nco(R) {2+ ^ + < la 2 n cstl b c^ + ^I~\ 



Proposition 5^3 follows directly from these choices plugged into the study of the 
cutoff case. □ 



5.2. Proof of Theorem |2.6[ Now we proved the immediate appearance of a lower 
bound only depending on the norm of the velocity we can spread it up to an expo- 



nential lower bound. As in Section 4.3, we thoroughly follow the proof of Theorem 
2.1 of [TH]. The proof in our case is exactly the same induction, starting from Propo- 
sition [573 and therefore we just describe briefly here how to construct the expected 
exponential lower bound. For more details we refer the reader to [10] . Section 4. 



As in Section 4.3, we spread the initial lower bound (Proposition 5.3 by induction 
where, at each step, we use the spreading property of the operator and fix e n 



small enough to obtain a strictly positive lower bound, see (5.8). 



There is, however, a subtlety in the non-cutoff case that we have to deal with. 
Indeed, at each step of the induction we choose an e n smaller and smaller, but at 
each step n the action of the operator — (Q~ + Q 2 ) behaves like the following (see 



(5.5)) 



exp 



-C f [m b Mco + n b co ) (t« - tW)(v) 



By (5.6) — (5.7), as e n tends to we have that n b co goes to +00 and therefore the 

1 1 ' 1 £n 

action of — (Q~ + Q 2 £ ) seems to decrease the lower bound to exponentially fast at 
each step. The idea is thus to find a time step A n = £„ — tn at each step to be 
sufficiently small so that it counterbalances the effect of n b co. 



More precisely, by starting from Proposition 5.3 as the first step of our induction, 
taking 

\k=0 / \fe=0 / 

in (5.5) and fixing e n by ( |5.8[ ) we can prove the following induction property 



Proposition 5.4. Let f be the mild solution of the Boltzmann equation described 
in Theorem \2.Q\ 

For all r in (0, T) and any sequence (A n ) ngN such that J2 n ^ A„ = 1, 



e N,Vt e 



J2 A 

vfc=0 



k]T,T 



,W(x,v) eUx R d ,f(t,x,v) ^ a n (T)l B (o,r n ){v), 



with the induction formulae 



a n+1 = est A n+ iexp 



[C/a^-^ 2 - 1 ]-* 



2 -Y+NtN/2+1 
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-cstlog^a^r^-W 2 - 1 ] ( £ r t 



fl 2 7+AVJV/2+1 
n n >n 



v is in (0, 2), 
a n+ i = est A n+ iexp 
if u — and 

r n+1 = V2r n (l - £ n ). 

(Cn)nGN any sequence in (0, 1) and r = ry, ao(r) and C/ depend only on r, Ef, 
E), W f (and L p f if $ satisfies §L4$ with 7 < 0). 



We emphasize here that the induction formulae are obtained thanks to the use of 



equivalences (5.6) and (5.7) inside the exponential term 

-C f (m bNCO +n hCO ) (4 X) -4 2 ) )<H)T f+ ^ -C/ (m^co+n b co) (Ek^„ + i A fc )(i^ + 

see step 2 of proof of Proposition 4.2, Section 4 in [TO] . 

As we obtain exactly the same induction formulae as in [TU], the asymptotic 
behaviour of the coefficients a n is the same. Thus, by choosing an appropriate 
sequence (A n ) n£ N, as done in [JO], we can construct the expected exponential lower 
bound as in the proof of Theorem 1.2, Section 4 in [TOj and make it independent of 



time in the same way we did at the end of Section 4.3 



Appendix A. Geometric study of grazing collisions 



This appendix is dedicated to the proof of Proposition 4.8 



We first recall that for x in U and v in M. d we defined, in Section |3| t min (x,v) 
to be the time of the first proper rebound when we start from x with a velocity 
—v. This means that t m i n (x,v) do not take into account when a ball rolls on the 
boundary and it implies that one cannot hope from continuity of the function t m i n 
as changing very slightly the velocity may lead to a proper rebound instead of a 
rolling movement. 

This being said we define a time of collision against the boundary which will not 
take into account the possibility of rolling along the boundary of hi. This will not 
be too restrictive as we are considering now a C 2 convex domain and therefore a 
trajectory that stays on the boundary will only reach a stopping point and this only 
happens on a set of measure zero in the phase space (see Section [3]). Therefore we 
define for x in U and v in IR d , the first forward contact with the boundary, t(x,v). 
It exists by the same arguments as for t min . Notice that if x is on dlA then for all 
v 7^ we have that t(x,v) = if and only if n(x) ■ v ^ 0, n(x) being the outward 
normal to dU at the point x. 

We have the following Lemma dealing with the continuity of the outward normal 
to dlA at the first forward contact point which will be of great interest for proving 



the crucial Proposition 4.8 



Lemma A.l. Let hi be an open convex bounded C 1 domain in S 
Then t : (x, v) — > t(x, v) is continuous from U x (W 1 — {0}) to 



Del 



) + 
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Proof of Lemma A.l Let suppose that t is not continuous at (xo, vq) in Ux (W 1 — {0}) . 
Then 

Be > 0,ViV ^ l,3(x N ,v N ), I V° „\\^] /at and \t(x ,v ) - t(x N ,v N )\ > e. 

I \\ v ~ V N\\ ^ l/iv 

If we denote by lu the longest diameter of U (it exists since the domain is 
bounded) we obviously have that for all N, ^ £(xjv,t>Ar) ^ Wll^ivll- Thus, 
(t(xjy, v N )) NeN is a bounded sequence of R and we can extract a converging subse- 
quence (t(ar^(jv), v<ft at))) such that T = lim t(x^ N ),v^ N )). 



By construction (see Section[3| we have that for all N in N, x<f>(N)+t( x <i>(N), v</>(N)) v (j>(N) 
belongs to dli which is closed. Moreover, this sequence converges to xq + Tvq which 
therefore is on dU. 

Finally we have that |t(xo,fo) — T\ ^ s and that U is convex. Therefore the 
segment [x ,x + max(t(x , v ), T)v ] stays in U and intersect the boundary at at 
least two distinct points. By convexity of the domain, this implies that the extreme 
points of the latter segment have to be on the boundary which means that xq belongs 
to dU which is not. 

Therefore, t is continuous inWx (M. d — {0}). By construction of t we have its 
continuity at the boundary (as n(x).v ^ means we came from inside the domain 
to reach that point and therefore the contact time tends to 0). □ 



Now we can prove the geometric Proposition 4.8 



Proof of Proposition |4~8} Consider e > and < v m < Vm- 

The first step is to understand that if a whole trajectory stays close to the bound- 
ary, then the angle made by the velocity with the normal at the point of collision is 
close to 7r/2 for dimension d = 2. The same behaviour in higher dimensions is de- 
scribed by the scalar product of the direction of the trajectory and the normal being 
close to zero. One has to remember that controlling ||Kj(x,f) — v\\ is the same than 
controlling scalar products with the boundary at each collision point (see definition 
of V s (x,v) in Section [3]). 

Take x on dU and p in N*. We define 

T p (x) = {\n(x) ■ v\ : v E § d_1 s.t. n(x) ■ v < and Vs G [0, t(x,v)], x + sv £ Ui/ P \ , 



with Uu p being defined by (4.6) 



r p (x) gives us the values of scalar products between a normal on the boundary 
and all the directions that create a characteristic trajectory which stays at a distance 
less than 1/p from the boundary in between two distinct rebounds (see Fig|5]). This 
is exactly what we would like to control uniformly on the boundary. 

We can remark that T p (x) is not empty because U and, thus, U\/ p are convex 
and by the geometric theorem of Hahn-Banach we can separate U\/ p and a disjoint 
convex containing x. It is also straightforward, a mere Cauchy-Schwartz inequality, 
that T p (x) is bounded from above by 1. Therefore we can define, for all p in N*, 
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hp: dU — ► R + 

x i — > sup T p (x). 

We are going to prove that (h p ) p( z^* satisfies the following properties: this is a 
decreasing sequence, h p is continuous for all p ^ 1 and for all x in dU {h p (x)) p£ ^* is 
convergent to 0. 

The fact that (h p ) is decreasing is obvious. 

In order to prove the continuity of h p we take an x on the boundary and v in S d_1 
such that \n(x) ■ v\ is in T p (x). We have that for all s in [0, t(x, v)] 

Vs 6 [0,t(x,v)\, d(x + sv,dU) < 1/p. 

The distance to the boundary is a continuous function and [0, t(x, v)] is compact so 
there exists s(x, v) in the latter intervalle such that d{x + s(x, v)v, dU) is maximum. 
Because U is convex we have that U\/ p is convex and therefore 

/ d(x + s(x,v )v,Ui/ v )\ 
Vs 6 [0, t{x,v)\, Blx + sv,^ v 2 ; l,p> J nZY Vp = 0. 

Then for all x' on the boundary such that \\x — x'\\ ^ d(x + s(x,v)v,Uy p )/2 we 



have that for all s in [0, t(x', v)], x' + sv is not in U\i v . Lemma A.l gives us that if x' 
is close to x then t(x',v) > and thus t> is not tangential at x either. Moreover U 
is C 2 so the outward normal to the boundary is continuous and therefore for x' even 
closer to x we have that v is such that \n(x') ■ v\ is also in T p (x'). To conclude, we 
just have to notice that the scalar product is continuous and therefore for all rj > 
we obtain, for x' even closer to x we have 



-v 



\n(x') ■ v\ — \n(x') ■ v\ 



The same arguments with the same constants (since our continuous functions act 
on compact sets and therefore are uniformly continuous) if x' is closed to x then 
taking \n(x') ■ v\ in Ti/ p (x') we have \n(x) ■ v\ in Ti/ P (x) and the same inequality as 
above. This gives us the continuity of h p at x by considering supremums in each 
inequality. 

Finally, it remains to show that for x on the boundary we have that h p (x) tends 
to as p tends to +oo. 

One can notice that the vector —n(x) is the maximum possible in T p (x) and is exactly 
the direction of the diametre in hi passing by x. Hence, simple convexity arguments 
leads to the fact that if all the segments of the form [x, x — t(x, —n(x))n(x)] intersect 
U\jp then we have that for all x on the boundary, it exists v p (x) in such that 
n(x) ■ v p (x) = —h p (x). Moreover, the segment [x,x + t(x,v p (x))v p (x)] is tangent 
to Ui/ P and we denote by x p its first contact point (see Figj5|. The convexity of U 
and Ui/ P shows us that x p is closer and closer from x and from the boundary (U is 
convex) and therefore v p (x) tends to a tangent vector of the boundary at x. This 
shows that 

lim h p (x) = 
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in the case where all the segments of the form [x,x — t(x, —n(x))n(x)} intersect U\i v . 

We now come to the case where the segments of the form [a;, x — t(x, —n(x))n(x)] 
do not all intersect tl\/ p . If for all p, this segment does not intersect U\/ p this implies 
by convexity of U that [x, x — t(x, —n(x))n(x)] is included in dU. But then —n(x) 
is not only a normal vector to the boundary at x but also a tangential one at x. 
Geometrically this means that x is a corner of dU and n(x) is ill-defined. This is 
impossible for U being C 2 . Hence, for all x on the boundary, it exists p(x) such 
that the segment at x intersect U p r x y However, U is C 2 and we also have Lemma 



A.l Those two facts implies that p(x) is continuous on dU which is compact and 
therefore p(x) reaches a maximum. Let us call this maximum P. For all p ^ P, 
all the segments of the form [x,x — t(x, — n(x))n(x)}, x in dU, intersect Up and we 
conclude thanks to the previous case. 

Thanks to these three properties we are able to use Dini's theorem, since dU is 
compact. We therefore obtain that (h p ) p ^* converges uniformly to 0. By taking 
p £ big enough we have that for a segment of a characteristic trajectory joining two 
points on the boundary to be outside U Ve we must have Y Ve ^ e for all x on the 
boundary (see Fig|5j). 



v 



Now we turn to more general trajectories. We take x in dU and v such that 
i ^ v ^ v m and we suppose that for a given t > 



Vse[0,*], X.(x,v) i U 1/P(e/2) . 

We are about to find a uniformly small time such that trajectories having at 
least two collisions against the boundary do not undergo an important evolution 
of velocity. This will be achieved thanks to the facts that \\v\\ ^ % and that the 
maximum of the scalar product is attained at a critical vector and this latter is the 
only one to be controlled. 

The first property to notice is that if (X s (x,v)) se [ 0jt j has at least two rebounds 
against the boundary, then at each of them the scalar product between the incoming 
velocity and the outward normal is less than 2e. Secondly, U is C 2 and therefore n(x) 
is uniformly continuous on the boundary. Thus, the specular reflection operator TZ X 
is uniformly continuous on dU x £>(0, Vm)'- 

3a>x > 0, Vx, x' G dlA s.t. ||a; — x'\\ < a x , \\TZ X — 1Z X >\\ ^ e/2. 



Finally, Lemma A.l tells us that (x,v) h- >■ x + t(x,v)v is unformly continuous on 
the compact dU x B(0, vm)- Gathering those last three properties and the fact that 
the lenght of the critical characteristic trajectory is greater than l/(p E / 2 ) raises that 
for 

{ a x 1 \ 

t ^ t e (v M ) = max — , 

\V M Ps/2V M J 

if (X s (x, v)) se [o t t] collides at least twice with the boundary then 

Vs G [0,t], ||v - V s (x,v)\\ ^ e. 
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Roughly speaking we do not allowed the velocities near the critical direction to 
bounce against the wall and for the grazing ones we run them for a short time, 
preventing them from escaping a small neighbourhood where the collisions behave 
almost the same everywhere (see Fig{5]). 

To conclude our proof, it only remains to find I ^ l/p e /2 that prevents trajectories 
staying in U\ but go through only one rebound with a scalar product greater than 
e/2 from happening. This is easily achieve by taking / small enough such that 
not a single trajectory with a scalar product greater than e/2 can stay inside IA\ 
during a time r. Indeed, one part of these trajectories will overcome a straight 
line of lenght at least v m r/2 and making a scalar product greater than e/2. The 
distance from the boundary of the extremal point of these straight lines is therefore, 
by convexity, uniformly bounded from below (e.g. in dimension 2 it is bounded by 
v m re/4:. Taking l e (v m ,r) being the minimum between this lower bound and l/p £ /2 
gives us the required distance from the boundary. 



X 




FIGURE 5. Control on grazing trajectories 

□ 
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